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ABSTRACT 

Micron size extrasolar dust particles have been convincingly detected by satel- 
lites. Larger extrasolar meteoroids (5-35 microns) have most likely been detected 
by ground based radar at Arecibo and New Zealand. We present estimates of the 
minimum detectable particle sizes and collecting areas for both radar systems. 
We show that particles larger than ~ lO/xm can propagate for tens of parsecs 
through the interstellar medium, opening up the possibility that ground based 
radar systems can detect AGB stars, young stellar objects such as T Tauri stars, 
and debris disks around Vega-like stars. We provide analytical and numerical 
estimates of the ejection velocity in the case of a debris disk interacting with a 
Jupiter mass planet. We provide rough estimates of the flux of large microm- 
eteoroids from all three classes of sources. Current radar systems are unlikely 
to detect significant numbers of meteors from debris disks such as Beta Pictoris. 
However, we suggest improvements to radar systems that should allow for the 
detection of multiple examples of all three classes. 

Subject headings: 



1. INTRODUCTION 

In astronomy, as in everyday life, most of our information comes to us in the form 
of electromagnetic radiation. Some astronomical systems also emit solid particles, which 
could yield valuable information if detected. The dust detectors on the Ulysses and Galileo 
spacecraft have yielded preliminary information on the flux of interstellar grains passing 
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through the solar system (Prisch et al. 1999; Landgraf et al. 2000). Due to the hmited 
area of the collecting surfaces, of order 200 cm^, these authors were unable to report fluxes 
for grains with masses > 10~^° g. These low-mass grains cannot be traced to their point of 
origin, since interstellar influences (e.g., gas drag, radiation pressure, magnetic fields) rapidly 
adjust their velocities. 

Particle fluxes probably decline dramatically as the grain mass increases, necessitating 
dust detectors with much larger collecting areas. Meteor-tracking radar facilities can be 
used to detect grains originating beyond the solar system, if the grain's initial velocity 
can be inferred with sufficient accuracy. The effective collecting area in this case could be 



Recently, Baggaley and coworkers reported the detection of extra-solar grains using the 
Advanced Meteor Orbit Radar, or AMOR (Baggaley et al. 1994; Baggaley 2000). Baggaley 
(2000) finds a "discrete source" with an angular diameter of ~ 30°, which he tentatively 
associates with /3 Pictoris. Inspection of his Fig. (2) suggests the presence of a broad, 
band-like feature as well. Baggaley does not provide any estimate of the particle fluxes. 

Our goal here is to identify likely sources of extra-solar grains and estimate expected 
fluxes at Earth, as a function of the grain size. We consider three types of objects that could 
potentially yield signiflcant fluxes: young main sequence stars, asymptotic giant branch stars 
(AGB stars), and young stellar objects (YSOs). 

For a given source of interest, the flux at Earth depends on the following three factors: 
1. The "dust luminosity" of the source, i.e., the rate at which grains are emitted. We also 
need to know whether or not the emission is isotropic. If not, then we need to know the 
orientation of the source. 2. The distance between the source (at the time when the grains 
were ejected) and the Sun (now). 3. The probability that the grain survives the trip and 
is not deflected on its way. Although deflected grains can be detected at Earth and will 
contribute to the general dust background, they do not reveal their source. 

Deflne the "speciflc dust luminosity" L^^o(t, t'ej, a) by 



where t is the age of the source, L(t) is the dust luminosity, v^j is the speed at which grains 
are ejected, and a is the grain radius.^ In order to calculate the dust flux at Earth, it is most 
convenient to consider the reference frame in which the source is stationary. We assume 



> lO^kml 




(1) 



•^Throughout this paper we make the simphfying assumption that the grains are spheres. 
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that, in this frame, large grains simply travel radially outward; i.e., we ignore the Galactic 
potential and gravitational interactions with individual stars, as well as any other influences 
that could deflect the grains (see §3.2). The number density of grains at distance d from the 
source is given by n{t, d) — J dv^j dariy^ait, v^j, a, d), where 

/, j\ Lvj,(^l\,y(l.f^,j f '^^/'•'cj) r r / rl\ • ( 0\ 

n'v,a(J'i '^eji O'l d) Aird^V /beam /survive l^^ej ) Ct) J > \-^) 

tej is the age of the source when the grains are ejected. The factor /beam accounts for 
anisotropic emission from the source (we assume that it is independent of Vej and a) and 
/survive ('f^ej) d) is the fraction of the grains that survive the trip out to distance d without 
being destroyed or significantly deflected. The particle flux at Earth at the present time is 
given by F{t) = / dvej da Fy^a{t, fej, a, d), where 

P /. i\ _ (. i\ _ Lv^a{Vej, a, te^ = t — d/ Vej) Vd,Q r r / i\ 

^Vjay^j "^ejj O, U) — n^Vjay^j ^ejj UjVd^Q — 47r(^2 y /beam /survive v'^^ej; O, U) 

(3) 

The velocity of the dust grain with respect to the Sun is given by 

Vd,Q = V^,Q - Vef^^Q , (4) 

where v^,^q is the velocity of the source with respect to the Sun and f*,© is the unit vector 
pointing from the Sun to the source. 

In §2, we estimate the limiting particle flux that could potentially be detected with 
radar facilities. In §3, we consider the mechanisms that can prevent a grain from reaching us 
and estimate the (size-dependent) distance that a grain can travel and still reveal its source. 
We estimate observable dust fluxes from young main sequence stars, AGB stars, and YSOs 
in §§4, 5, and 6, respectively. We discuss our results in §7, and summarize our conclusions 
in §8. 



2. DETECTABLE FLUXES 

2.1. Satellite Measured Fluxes 

There have been several claimed detections of extrasolar meteors, both from satellite 
detectors and ground based radar. The most convincing detections are those of the Ulysses 
and Gahleo satellites. Prisch et al. (1999) and Landgraf et al. (2000) show that the flux 
of m 6 X 10~^^g extrasolar particles is given by mfm ~ 10"^ cm~^ s~^. Here fmdm is 
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the flux of particles with masses between m and m + dm^ At the low mass end of the 
mass distribution (below m ~ 6 x 10~^^ g) the particles are subject to strong perturbations 
from the solar wind and its associated magnetic field, so that the measured flux is not 
representative of the flux of small particles outside the heliosphere. Above m ^ 6 x 10~^^ g, 
they find mfm oc m~^^^] there are roughly equal masses of particles in every logarithmic mass 
bin. 

For purposes of extrapolation we will occasionally assume that the cumulative flux of 
interstellar particles of mass m is 



for m > 6 X 10 g. 

We employ two other differential fluxes. The flrst is /„, where fada is the flux of particles 
with radii between a and a + da\ fa has units cm^^ s^"*^ cm^-*^. In many cases one assumes 
that fa follows a power law distribution, oc a^"' . This corresponds to oc m~", with 
q; = (7 + 2)/3. We also use the differential mass ffux, m^fm, with units gcm~^s~^. 

In the asteroid literature one encounters the Dohnanyi (1969) law, fm oc m~", with 
a = 11/6, corresponding to 7 = 7/2. The Dohnany law represents a singular, steady state 
solution to a set of equations describing a closed system of colliding bodies. Asteroids do 
not follow this relation (Ivezic et al. 2001). The same scaling appears in the interstellar dust 
hterature, in the MRN dust model (Mathis, Rumpl, & Nordsieck 1977). In these models the 
mass in large particles diverges in proportion to ^/a, or as m}^^. The observed scaling (eqn. 
5) is a 2.1, shghtly larger than the Dohnanyi value, so that the mass in large particles is 
finite. We note that particles in debris disks or in the outflows of AGB stars need not be 
distributed according to either law; in fact the findings below suggest that they are not well 
described by a power law with j — 7/2. 



Ground based radar at Arecibo Meisel, Janches, & Mathews (12) and in New Zealand 
(AMOR) Taylor, Baggaley, & Steel (18); Baggaley (1) have also reported detections of 
extrasolar meteors. When a meteoroid enters the Earth's atmosphere, air molecules ablate 



''The cited authors use the cumulative flux in a logarithmic mass interval, /iogiom£^(logTO) = fmdm. The 
two fluxes are related by mfm = /logm m/ In 10. 





2.2. 



Radar Fluxes 
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and ionize material from the meteoroid. The free electrons created by this ablation reflect 
radio waves, a fact exploited by radar aficionados. 

The size of individual meteoroids detected by radar systems can be inferred by three 
different methods; by the radar power reflected from a Presnel zone of the meteor trail 
(§2.2.3), by the power reflected from the meteor head (§2.2.4), or by the deceleration of the 
meteoroid as it ablates in the Earth's atmosphere (also in §2.2.4). 

The first two size estimates rely on the power of the reflected radar signal. For radar 
micrometeors the relevant reflection is coherent. The wavelengths employed are comparable 
to or larger than the initial width tq of the ionization trail, so the radar is sensitive to the 
electron line density q rather than the space density. The reflected power depends on the 
transmitted power, the distance between the radar and the meteoroid, the mass and velocity 
of the meteoroid, and the (uncertain) ionization efficiency. Here we quantify this relation. 
A good introduction to the subject of radar meteors is given by McKinley (14). 



2.2.1. Meteoroid properties 

The radii of radar meteoroids range from a ~ 3 /im to ~ 40 jim.. For purposes of 
illustration, we consider meteoroids with radius a = 10 /xm, density p = 3gcm^'^, and made 
of silicates such as Forsterite, Mg2Si04. Our default meteoroid has a mean molecular weight 
of 140 grams per mole, and a vaporization temperature of order 2000K. We assume that the 
energy per bond is ~ 2 eV. This is equivalent to a heat of ablation of C ~ 10^^ ergs per gram. 

The mass of our default meteoroid is m = 1.3 x 10^^ g. The mean atomic weight is 
H PS 20mp, where rrip = 1.67 x 10~^^g is the mass of a proton, while the total number of 
atoms in the meteoroid is = m/jj, pa 4 x 10^"^. The total binding energy of the meteoroid 
is 

A typical meteor is highly supersonic, with velocity w ~ 40 km/ s at the top of the 
atmosphere, so the kinetic energy of the meteoroid is 




much larger than the binding energy. 
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2.2.2. The ionization trail 

As the meteoroid passes through the atmosphere, it colhdes with air molecules and is 
ablated. Because the kinetic energy per meteoroid atom is in excess of 100 eV, while the 
ionization potential is of order lOeV (depending on the atomic species), collisions between 
ablated atoms and air molecules may ionize either or both particles. The number of ions 
produced by each meteoroid atom is denoted by (5 in the meteor literature. It is electrons 
from these ionized atoms that reflect the radar signal. 

If the mass loss rate of a meteoroid traveling at velocity v is dm/dt, the number of ions 
produced per centimeter along the flight path is 

g = (8) 

fiv dt 

The rate at which the meteoroid ablates depends on the ambient air density, which increases 
exponentially as the meteoroid descends through the atmosphere. In other words, the mass 
loss rate increases rapidly with time. The pressure scale height Hp of the atmosphere is 
roughly constant, since Hp = kT/diaird) and the air temperature T ^ 200 K. The line 
density q reaches a maximum just before the meteor vanishes, with most of the matter being 
deposited in the last scale height of its path. Taking —vdt — Hp, we find 

3 /3m 

where the factor of 3/2 comes from the more accurate calculation outlined in the appendix. 

The value of (3 depends strongly on velocity v, but the range of v for typical extrasolar 
meteors at Arecibo or AMOR is 20 — 60kms~^, with most near 40kms~^. Jones (1997) 
and Jones & Halliday (2001) examine (3 using a combination of laboratory experiments and 
observational data. For iron atoms they find /3 ^ 0.6(f/40kms^^)^'^^ (with v < 60kms~^). 
For more realistic compositions they expect to be a factor of 5 smaller. For a meteoroid 
with mean atomic weight = 20mp (as for Forsterite), 

, 6 X 10^ f ( P ) ( em-. (10) 

The height at which the meteoroid ablates can be estimated by finding the height at 
which the binding energy of the meteoroid decreases most rapidly. The rate of change of the 
binding energy is 

= -A . Tra^ . (11) 
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where pa is the mass density of the atmosphere at the height of the meteoroid and A is 
the fraction of the kinetic energy flux that goes toward ablating the meteor. Taking the 
derivative of this expression with respect to time and setting the result equal to zero yields 
the relation 

2 dm 1 dpa Sdv , , 

— _| (12 ) 

3m dt Pa dt V dt 
We show in the appendix that the last term is typically smaller than either of the other two, 
or at best comparable, so we ignore it. Using equation (11), and setting dpa/dt — pav/Hp, 
we find that the meteoroid ablates when the atmospheric density is 

Using typical values, we find 

1 3 10-9 ^ P yiO-^U « U C W 40kms-^ y /6km\ 

(14) 

Using an MSIS model atmosphere^ we find that this density occurs at a height of 96 km, 
comparable to the observed height at both AMOR and Arecibo. 

On the other hand, radar data from Jodrell Bank and from Ottawa show that radar 
meteoroids detected by those systems also ablate at heights near 95 km (McKinley 1961), de- 
spite the fact that they have line densities q ~ 10^^ cm"-*^, six orders of magnitude larger than 
the Arecibo or AMOR meteoroids traveling at the same velocity. These larger meteoroids 
are still much smaller than the mean free path, so that aerodynamic effects seem unlikely to 
explain the difTcrcncc; nevertheless, something must differ between the two classes of objects, 
else the more massive Jodrell Bank meteoroids would be ablated five scale heights below the 
less massive meteoroids discussed here. One possibility that has been suggested is that the 
more massive meteoroids (those seen at Jodrell Bank and Ottawa) fragment into smaller 
pieces, which are then rapidly ablated. The observed heights of the Arecibo and AMOR 
events suggest that for micrometeoroids that do not fragment A 0.1 and ( ^ 10^^ cm^ s~^. 

We stated above that the initial width ro of the ionization trail was comparable to or 
smaller than the radar wavelength. We estimate ro as follows. The cross-section of an atom 
is roughly a = ird^, with d ^ 2 x 10~^cm; at 95km the MSIS atmosphere has T = 171K, 
and a mean free path of / ~ 37 cm. However, the typical meteoroid atom has an atomic 
mass about twice that of a nitrogen atom. In a collision between a meteoroid atom and 



^Available at http://nssdc.gsfc.nasa.gov/space/niodel/models/msis.html 
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a nitrogen molecule, the momentum will not be shared equally between the two N atoms; 
the binding energy of the N2 is much less than the kinetic energy of the meteoroid atom. 
Neglecting the momentum carried away by the slower of the two rebounding N atoms, we 
find that the final velocity of the meteoroid atom is about 2/3 of its initial velocity for a head 
on collision. After n such collisions, the velocity is f„ = (2/3)"fo. Setting f„ equal to twice 
the thermal velocity at 95 km, we find that the meteoroid atom is thermalized after n ~ 10 
collisions. Thus Tq ~ y/nl^ or about 100cm. A strict upper limit to the thermalization time 
is Tth — l/v-i about 10/xs, much shorter than the time between radar pulses. Using different 
arguments, (Manning 1958) and (Bronshten 1983) obtain similar results. After this initial 
rapid diffusion, ordinary thermal diffusion sets in and the trail radius slowly increases. 

The initial electron density is 

- - - - X (n) (li;)' (^) (3^) 

The plasma frequency ^ 0.39 MHz is well below the frequency employed in the radar 
systems, so the radar beam will penetrate through the trail, ensuring that all the electrons 
will reflect. The average distance between electrons is n~^/^ ~ 0.1 cm, much less than the 
wavelength of either AMOR or Arecibo, so the reflected emission will be coherent in both 
systems. 

The thermal diffusion mentioned above eventually causes the trail radius to exceed the 
wavelength of the radar. This diminishes the visibility of the trail, particularly for trails 
formed at great heights, or for short wavelength radar. A recent discussion of this 'height 
ceiling' can be found in Campbell-Brown & Jones (2003). 

2.2.3. Minimum detectable meteoroid size at AMOR 

We proceed to calculate the power collected by the radar receiver at AMOR. The relevant 
properties of the AMOR radar (as well as those for the Arecibo setup) are listed in Table 
1. An electron a distance R from the radar sees a flux PtGt/ -^t^R^ , assuming the radar 
transmits a power Pt and has an antenna gain Gt- The gain is defined as 47r steradians (the 
solid angle into which an isotropic emitter radiates) divided by the solid angle of the radar 
beam; the latter is roughly (X/D)^, where D is the radius or length of the antenna. Hence the 
effective area of the antenna is related to the gain by A^fj = CxX'^/^tt. The electron scatters 
the incident wave, emitting a power {^/2)PtGtcft / ^'kR^ back toward the radar; the factor 
3/2 arises because Thompson scattering is that of a Hertzian dipole rather than isotropic. 
The receiving antenna (which in the case of Arecibo is the same as the transmitting antenna) 
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captures a fraction Aeff/AnR'^ of this scattered power. The power received at the antenna 
from a single electron is 



i?4 



(16) 



1287r3 

We have allowed for the possibility that the gain for the receiver and transmitter are different, 
as is indeed the case at AMOR. 



The amplitude of the electric field at the receiver is 
where r is the impedance of the receiver. 



(17) 



The instantaneous signal amplitude received from a meteor trail is found by summing 
over all the electrons along the track: 



A{t) 



(2rAPe)^^^ q{s) sm{ujt - 2kR)ds, 



(18) 



SI 



where lj is the angular frequency of the radar, and k = 27t/X. The total distance from 
the radar to the trail and back is 2R, whence the factor 2 in the argument of the sine. 
We measure s along the trail starting from the point X on the trail at which q reaches its 
maximum. The distance between the radar and point X is denoted by Rq, while the angle 
between the hne of sight and the meteor trail at X is denoted by 7r/2 + ^. 



Near X, 



R^ Ro 



S 
Rq 



(19) 



where we assume 9 << 1. Define % = a;t — 2kRo, and z = 2s/ \/Rq\. Then 

v2 



A{t) = V2rAPeg 



sm 




Rq TTZ- 

^_2n\l-ez- — 



dz 



(20) 



When 6 « \J\/ Rq the integral reduces to a Frcsnel integral, which is an oscillating quantity 
of order unity. The peak power at AMOR is then 



T" 



2567r3 VP, 



q ar- 



(21) 



When 9 > a/A/Pq, the integral cuts off at 2; < 1/4. In other words, the AMOR detector is 
sensitive only to micrometeor trails within an angle 9 ~ -^/A/Po of the perpendicular to the 
hne of sight. 
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The minimum detectable line density for AMOR is 

1 /Pfl\"V 25671^ \"Vfio"''" 



The minimum detectable power at AMOR is P„ = 1.6 x 10^^"^ Watts, the transmitted 
power is Pt = lOOkW, while the gains are Gt = 430 and Gr = 130. The typical range is 
Rq pa 200 km, leading to 

q„nn - 7 X 10^ cm-\ (23) 
Using equation (9), the minimum detectable mass is 

rrimin ~ 2 X 10-' g, (24) 

corresponding to a minimum radius of 

o -3\ 1/3 /n 1 \ 1/3 / D \ 2/3 

3gcm'^\ /0.1\' / /to ^ 



We are now in a position to estimate the collecting area A^oii of the AMOR detector 
(as a meteor detector, not as a radar). The collecting area is the product of the sensitivity 
weighted geometric area Aq imaged by the radar and the fraction of solid angle ^y X/Rq to 
which it is sensitive; 

Acou^J^Ag. (26) 
We define the sensitivity weighted geometric area as 

Here F(m) is the flux of particles of the minimum detectable mass m — mmin{R)] the latter 
is a function of i?, the distance between the radar and the meteor trail. We assume that 
T^min oc qmin oc R^^"^, whcrc the last scaling comes from equation (22), and that F{m) ~ m"^ 
as an analytically simpler version of (5). With these assumptions, 

Ag ^ 2Rl,^A<P [VRmax/Rmin ' l] (28) 

The AMOR beam has a width of A0 3°, and extends from 12° to 30° above the horizon. 
We further assume that all the meteors are detected at heights of ~ 90 km. This implies 
that the lower and upper limits to the integration are Rmin = 165 km and Rmax = 430 km. 
We find Ag ~ 10^^ cm. The collecting area of the radar system is then 



^col 



8km^ (29) 
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Prom Fig. lb in Baggaley (2000), we estimate that in 4 years of continuous operation 
AMOR detected ~ 10^ extrasolar meteors, so we estimate the flux as 

mfm (m = 2 X 10"^ g) fti 4 x 10"^^ cm^^ s'K (30) 

From the same paper, we estimate that the number of extrasolar meteors detected from 
Baggaley 's "point source" to be about about 200, leading to a flux 50 times smaller. Both 
fluxes are plotted in Fig. (1), along with flux estimates from Ulysses and Galileo, and from 
Arecibo. 

Fig. (2) shows the same data as a differential mass flux. In this plot the bulk of the 
mass flux is represented by the highest point, making it easy to see that most of the mass of 
interstellar meteoroids reaching the inner solar system comes in the form of ~ 0.3/im objects. 
On the same figure we have plotted the size distribution found by Kim & Martin (1995), 
normalized to a hydrogen number density of 0.1 and assuming that half the metals are in 
the form of dust grains. 

2.2.4- Minimum detectable meteoroid size at Arecibo 

The minimum detectable meteoroid size for Arecibo can be calculated in a similar 
manner, starting from equations (16) and (17). Arecibo uses a much shorter wavelength 
than AMOR, 70 cm rather than 1145 cm. Arecibo also has a much narrower beam than 
AMOR, ~ 500 m versus ~ 10 km. It is thus insensitive to trails traveling across the beam, 
in contrast to AMOR. Due to the fact that trails observed perpendicular to the line of sight 
have long coherence lengths, AMOR is sensitive to fairly small meteoroids despite the rather 
modest power and antenna gains employed. Arecibo, by contrast, is sensitive primarily to 
vertical trails, which have short coherence lengths (a quarter of a wavelength, less than 
20 cm). This is compensated for by the much larger gains, and by the larger transmitted 
power of the Arecibo radar. 

The short wavelength employed makes the Arecibo radar sensitive to diffusion of elec- 
trons across the meteoroid trail. The Arecibo radar beam is roughly vertical, with a half 
opening angle of 1/6°, so at a height of 100 km it is about 500 meters wide. Typical Arecibo 
meteor trails are a few kilometers long, so they are within a tenth of a radian or so of the 
vertical. Thus the ion trail left by an Arecibo meteoroid diffuses primarily horizontally, since 
the density along the trail varies on kilometer scales. As shown above, the initial trail width 
is of order 100 cm, but this width increases over time due to diffusion. The diffusion coeffi- 
cient ~ C5Z ~ 3 X lO^cm^s"^, where Cg ~ Sx lO^cms""^ is the sound speed, and Z 10 cm 
is the mean free path of air molecules. In the horizontal direction the Fresnel length is 
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Lp = VRqX 3 X 10^ cm, much larger than the initial beam width tq ~ 100 cm. The time 
to diffuse horizontally over a Fresnel zone is L^/D k, 3000 s, so diffusion is irrelevant for a 
vertical trail. However, for a trail 0.1 radians off the vertical, the time to diffuse A/4 in the 
vertical direction is about 13 ms. This will limit the time the beam remains visible to the 
radar at a given height to about 13 ms for a typical meteor. This is consistent with the trail 
lengths observed by Mathews et al. (1997). 

We have already noted that the meteor trail is much narrower than the Fresnel length 
of the Arecibo radar; then all the electrons in a patch of length of order A/2 along the trail 
radiate coherently. Suppose that Nch is the number of coherently emitting electrons in such 
a patch. The amplitude of the electric field at the receiver from a single patch is 



A ^ V2rAPeA^c/.- (31) 

The vertical resolution of the Arecibo antenna is L pa 150 m, so that there are Np L/ \ k, 
200 patches radiating in a resolution length. The signal observed at Arecibo varies on a 
length scale much shorter than the atmospheric scale height, indicating that the mass loss 
rate fluctuates rapidly. We assume that these fluctuations occur on scales smaller than L, so 
that the electric held from each patch adds incoherently. The power received by the detector 
is 



\2 



^-0 



Since the coherence length along the (assumed vertical) trail is of order A, the number of 
coherently scattering electrons is simply related to the line density, 

N,h^q- A/2. (33) 

The expression for the power becomes 

For a given transmitted power, wavelength and antenna gain, this is smaller than the 
equivalent expression (21) by the ratio 



.10-^ (35) 



using values for Arecibo. Worse, the power scales as A^ (or as A^-^ for Arecibo), so the ~ 15 
times shorter wavelength employed by Arecibo reduces the received power by another factor 
of ~ 3000. However, these two factors are more than compensated for by the much larger 
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gains of the Arecibo antenna. The net effect is that despite its higher power and much larger 
gain, the Arecibo radar can detect meteoroids that are a factor of seven to ten smaller in 
radius than the AMOR setup, but no smaller. 

Define the signal to noise ratio SNR = P/Pn, where P„ is the noise power for the 
Arecibo receiver. Using this definition, and solving for 



Using the values in Table 1, we find q ~ 10^ cm ^. Using equation (9) we find 
The corresponding meteoroid radius is 

«-mly)lT^)l^)"-«-- ,33) 

A firm lower hmit for the size of the meteoroids seen at Arecibo is found by noting that 
the energy per meteoroid atom is ~ 100 eV. Since the ionization potential of either atoms or 
molecules is typically of order 10 eV, a single meteoroid atom can free at most 10 electrons. 
Taking p — 3g/ cm^ and /3 = 10, we find amin ~ l-3//m. 

Prom equation (34), and noting that GrX^ — const., we see that, as long as A is larger 
than the initial radius Tq of the trail, P oc Some Arecibo meteors are seen at both 430 
MHz and 50 MHz (Meisel, private communication). This suggests that the volume occupied 
by the reflecting electrons has linear dimensions smaller than X/2 ^ 35 cm. 

For example, it might be thought that the meteoroid itself captures enough electrons 
to scatter the radar beam. However, this is not possible. The number of charges required 
to reflect the beam is given by equation (32), except that in this case there is only one 
emitting region (the meteor) so the factor \/L/X should be dropped. We find Nch ~ 10^; 
the meteoroid is charged to a (negative) voltage of V = 4.8 x 10^(3/im/a) volts. Charging 
the meteoroid to such a high voltage leads to rapid electron emission by quantum tunneling, 
also known as field emission (Fowler & Nordheim 1928). The discharge current (in statamps 
per square centimeter) is given by 



2nh (f> + fJ> 
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where E ^ V/a is the electric field near the surface of the meteor, is the work function, 
and fjL is the Fermi energy. Both of the latter have values of order 5 eV. 

We make the optimistic assumption that the meteoroid gains an electron for every atom 
it collides with in the atmosphere, /j„ = en*!), also in statamps per square centimeter {n* is 
the number density of air molecules at the height where the meteoroid ablates, see equation 
14). Setting these two currents equal, we find the equilibrium field E f» 3.6 x 10^ Vcm~^ or a 
maximum voltage of V = 10''(a/3;um). The number of charges on the meteor cannot exceed 
Nch.max ~ 2 X 10'^(a/3yum)^, a factor of 30 too small to explain the Arccibo observations. 
The meteor head signal must be due to electrons in the atmosphere around or trailing the 
meteor. 

One may also estimate a dynamically, if one can measure the rate of deceleration of the 
meteoroid due to atmospheric drag; 



where F is the drag coefficient, and is the area of the meteor. The radar measures dv/dt, 
V, and the height of the meteor. Using a standard atmospheric model, the density pa can 
be determined, For micrometeoroids, which have a size less than the mean free path for 
collisions, F = 1, so the ratio m/Am can be determined directly from measured quantities. 
(Janches et al. 2000) refer to this ratio as the ballistic parameter, BP. We assume spherical 
meteoroids, so BP ^ Apa/2>. Then a simple estimate for a is 



(Janches et al. 2000) find a range of baUistic parameters, ranging down to 10 ^gcm ^ or 
slightly lower, leading to 



This is a factor 2 smaller than the radius estimate based on received radar power, or a factor 
of 8 smaller when considering the meteor mass. For the meteors with the smallest BP, the 
size drops to 1pm or less. It is very difficult to see how the power reflected from such a low 
mass meteoroid could be detected at Arecibo, suggesting that the micrometeoroids have a 
density closer to lgcm~^ than to 3gcm~^. In fact, if we set p = 1, we find a ~ S.5pm and 
a ~ 7.5pm from the power and BP estimates, respectively. 



m— = -TpaV^A, 



(40) 




(42) 



The fluxes observed at Arecibo are given by Meisel et al. (2002); they are also shown 
on Figure 1. 



- 15 - 



3. PROPAGATION OF LARGE GRAINS THROUGH THE ISM 

We have seen that AMOR can only detect grains with masses > 2 x 10^^ g, corresponding 
to a grain radius a ^ 25/xm for sihcate dust (mass density p ~ 3.5gcm~^). Similarly, 
Arecibo can detect particles with a > 6/im. Thus, we will restrict our attention to "large" 
(by interstellar dust standards) grains, with a > lO/^m. 

For such large grains in the interstellar medium (ISM), radiation pressure can be ignored. 
The force due to radiation pressure is Frad = 7ra^((5pr)AMrad, where (Qpr) is the radiation 
pressure efficiency factor averaged over the interstellar radiation field (ISRF) and Awrad is 
the energy density in the anisotropic component of the ISRF. Adopting the ISRF for the 
solar neighborhood (Mezger, Mathis, & Panagia 1982; Mathis, Mezger, & Panagia 1983) 
and a 10% anisotropy (Weingartner & Draine 2001b), (Qpr) ~ 1 for a > 10//m and Au^ad ~ 
8.64 X 10~^^ergcm~^. The time interval required for radiation pressure to change a silicate 
grain's velocity by 1 kms~^ (assuming the anisotropy direction remains the same as the grain 
moves through space) is thus given by At ~ 5 x 10*^ yr(a/10/im). For large grains in the 
diffuse ISM, the forces resulting from the asymmetric photon-stimulated ejection of electrons 
and adsorbed atoms are even less effective than the radiation pressure (Weingartner & Draine 
2001b), so these forces can be neglected as well. 

The following three influences can be significant for grains with a > lOyum: 1. the 
drag force, 2. grain destruction following impacts with interstellar grains or high-energy gas 
atoms, and 3. the magnetic force. 

We will consider the propagation of grains through two idealized phases of the ISM, 
the cold and warm neutral media (CNM and WNM, respectively). We also consider the 
Local Bubble (LB), i.e., the large volume of low-density, ionized gas that surrounds the Sun 
(see, e.g., Sfeir et al. 1999). Our adopted values for the gas temperature Tgas, H number 
density rin, and electron fraction Xe = JT-e/^H (^e is the electron number density) in these 
three environments are given in Table 2. In each case we adopt the ISRF for the solar 
neighborhood. 

The drag and magnetic forces both depend on the grain's electric charge, which is set 
by a balance between the accretion of electrons from the gas versus photoelectric emission 
and proton accretion. For large enough grains, the rates of each of these processes are 
proportional to the grain area; thus, the electric potential $ of large grains is independent of 
the size a. Using the charging algorithm from Weingartner & Draine (2001c), we find that 
silicate grains with a > lO^m charge to (0.15 V, 0.74 V, 0.94 V) in the (CNM, WNM, LB).^ 



^Here, and throughout this paper, we take silicate optical properties from Li & Draine 2001. 
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3.1. The Drag Force 

When a grain's motion is highly supersonic, the hydrodynamic drag (i.e., the drag due to 
direct impacts of gas atoms and ions) does not depend on the gas temperature. In this case, 
the grain's speed decreases by a factor of e once it has encountered its own mass in gas. The 
sound speed is 0.81 km s~^ (6.2kms~^) in the CNM (WNM), assuming that the He number 
density riHe = O.lnn. We expect grain speeds as low as ~ lOkms"^, so deviations from the 
highly supersonic limit can be expected in the WNM. In addition to the hydrodynamic drag, 
there is also the Coulomb drag, due to long-range electric interactions between the charged 
grains and gas-phase ions. The Coulomb drag contributes significantly to the total drag in 
the WNM for drift speeds < lOkms^^. Using the approximate drag expressions from Draine 
& Salpeter (1979), we find that the actual drag force (hydrodynamic plus Coulomb) in the 
WNM exceeds the hydrodynamic drag force calculated in the highly supersonic limit by only 
a factor of f=:i 2 at a drift speed of lOkms"^. For regions outside of the LB, we will adopt 
the highly supersonic limit and assume that drag limits the distance a grain can travel to 
its velocity e-folding distance D^iag, given by 

Near the Sun, the average H density near the Galactic midplane 1 cm ^ (Whittet 

1992). 

In the LB, the sound speed is lOOkms"^; thus, the highly subsonic limit may apply. 
In this case, the drag force is proportional to the grain speed. Using the Draine & Salpeter 
(1979) drag expressions, we find that the Coulomb drag is negligible and that 

where Vq is the grain's initial speed. Since the maximum extent of the LB is ~ 250 pc (Sfeir 
et al. 1999) and tih ~ 5 x 10~^ cm~^, drag in the LB is always insignificant, whether the 
subsonic or supersonic limit applies. 



3.2. The Magnetic Force 

Charged grains spiral around magnetic field lines with a gyroradius 
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Although the magnetic field strength B has not been measured in the Local Bubble, B ~ 
just outside the LB, with the random component dominating the ordered component (Heiles 
1998). Thus, the magnetic force can significantly defiect grains with a ~ lO/xm. Since 
rs oc a^, the importance of the magnetic defiection rapidly decreases with grain size. 



When a grain travels through the ISM, it is subjected to various agents of destruction, 
including: 1. Sputtering, in which gas-phase ions strike the grain and remove grain surface 
atoms, which then enter the gas. It is useful to distinguish between thermal sputtering (due 
to thermal motion of the gas atoms) and non-thermal sputtering (due to the motion of the 
grain with respect to the gas). 2. Shattering and vaporization in grain-grain collisions. 3. 
For ices, sublimation and photodesorption. 

Thermal and non-thermal sputtering have been extensively discussed by Ticlcns et 
al. (1994). Applying their analysis, we find that sputtering of grains with a > 10/xm can be 
neglected in all environments of interest. 

The physics of grain-grain collisions has been treated extensively by Tielens et al. (1994) 
and Jones et al. (1996). When the relative speed exceeds ^ 3kms^^ (see Table 1 in Jones et 
al. 1996), a crater forms on the larger grain (the target). A portion of the evacuated mass 
remains on the grain as a crater lip. The rest of the crater mass is removed, both as vapor and 
shattering fragments. The shattering fragments dominate the vapor, and the largest fragment 
is typically somewhat larger than the smaller impacting grain (the projectile). When the 
relative speed exceeds ~ lOOkms"^ (depending on the target and projectile materials, see 
Table 1 in Jones et al. 1996), the larger grain is entirely disrupted, and the largest shattering 
fragment can be a substantial fraction of the target size. For our purposes, the cratering 
regime is generally applicable. 

Jones et al. (1996) give simple approximations for the (material-dependent) ratio of the 
crater mass (Mc) to the projectile mass (Mproj)- Interstellar dust is thought to be dominated 
by two populations: silicates and carbonaceous (e.g., graphite) grains. For a head-on collision 
between a silicate target and cither a silicate or graphite projectile, the Jones et al. (1996) 
results can be approximated by 



where v is the grain-grain relative speed. This empirical approximation reproduces detailed 
calculations, using equation (1) and Table 1 from Jones et al. (1996), to within 10% when 



3.3. 



Grain Destruction 




(46) 
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5 < V < lOOkms-i. The crater mass is approximately 4.2 (2.2) times bigger for a graphite 
(ice) target. We will assume that equation (46) applies to all collisions and that the entire 
crater mass is ejected. This is a conservative assumption, since the crater mass should 
actually be smaller for oblique collisions, and a portion of the crater mass will form a lip. 
On the other hand, we underestimate destruction if target grains are made of less resilient 
material than silicate or if the grains are fluffy. 

Since the mass in dust is ^ 0.011 times the mass in H (in all gas-phase forms) in the 
ISM (see, e.g., Wcingartner & Draine 2001a), the distance a grain can travel before it loses 
half its mass is given by 




(47) 

yielding Ddest ~ 860 pc for v — lOkms ^ (and with other canonical parameter values as in 
eq. 47). 



3.4. The Survival Probability 

In Figure 3, we plot D^rag, t^b, and D^est versus v for four grain sizes, covering the range 
of interest. Due to the expected dramatic decrease in flux with a (e.g., see eq. 92 below), we 
are not interested in a > 100//m. We adopt p = 3.5 gcm"^, riu = 1 cm"^, $ = 0.5 V, and 
B = 5 fiG, and assume supersonic grain speeds, for these plots. For all of the considered 
grain sizes, magnetic deflection dominates the other two processes when v = lOkms^^, 
while destruction dominates when v > a few to several xlOkms"^. The drag force is never 
dominant. The detailed analysis of the dynamics of a charged grain in a region with both 
ordered and random magnetic field components is beyond the scope of this paper. Thus, we 
will simply assume that 

J. _ r 1 , d < min(rB, -Ddest) 
1 ,d> min(rB, D^est) 

The grain velocity used in evaluating tb and -Ddest should be taken with respect to the 
ambient ISM. However, for simplicity, we will use the velocity with respect to the Sun. 

For Local Bubble conditions, magnetic deflection dominates destruction for grains with 
10 < a < 100/xm when v < lOOkms"^. For example, for a = 10/ym, Ddest = t'b when 
V — 311kms~^ and Tb = 281 pc (which exceeds the maximum extent of the LB). For 
a — 100/xm, Ddest — fB when v = 143kms~^ and rs = 13kpc. 



(48) 
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4. DUST FROM YOUNG MAIN SEQUENCE STARS 

IRAS observations of several young main sequence stars revealed emission at 60 and 
100/im substantially in excess of the emission from the star's photosphere (e.g., Aumann 
et al. 1984, Gillett 1986). This "Vega phenomenon" (named for the first example to be 
observed) was attributed to circumstellar dust, which absorbs the star's optical radiation and 
re-emits in the infrared. Shortly following this discovery. Smith & Terrile (1984) observed 
the optical hght scattered by the dust around P Pictoris and found that the grains lie in an 
edge-on disk. More recently, disks have been imaged (in the infrared and sub-mm) around 
a handful of other stars (e.g., Holland et al. 1998, Schneider et al. 1999). For reviews of the 
Vega phenomenon and circumstellar dust disks, see Backman & Paresce (1993), Lagrange, 
Backman, & Artymowicz (2000), and Zuckerman (2001). It is not yet clear what fraction 
of Vega-like stars posses dust disks, but if these stars also posses planets, then dynamical 
interactions can eject the grains from the system. (Radiation pressure also removes small 
grains, but these are too small to be traced back to their source.) 

The most massive disks, which orbit the youngest stars, have optical depths in the near 
IR of order 10^^ to 10^^. The bulk of this optical depth is contributed by grains with size 
oo ~ A/27r, where A is the wavelength, typically 1 — lO/xm (see below). 

4.1. Gravitational ejection of small particles 

Jupiter mass planets are seen around 5 — 10% of nearby solar type stars in radial velocity 
surveys. These surveys are not sensitive to planets in orbits larger than ~ 5 AU, such as 
Jupiter, so the fraction of solar type stars with such planets is likely to be substantially 
higher. Observations of debris disks also hint at the presence of planets (SchoU, Roques & 
Sicardy 1993; Wilner et al. 2002). This suggests that gravitational interactions between a 
massive planet and dust particles in debris disks are a natural means for producing interstellar 
meteoroids. 

We estimate the ejection velocity of large (25/xm) dust grains interacting with a Jupiter 
mass planet. We neglect collisions with other dust grains, an assumption which we justify 
at the end of the calculation, as well as radiation pressure. We follow the derivation of Opik 
(1976); a good general introduction to the two body problem can be found in Murray and 
Dermott (1999). 

We assume that the planet, of mass Mp, orbits the star, of mass M*, on a circular orbit 
of semimajor axis a^. A test mass (the dust particle) also orbits the star. The test mass 
is subject to the gravity of both the planet and the star, but the planet is assumed to be 
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immune to the gravity of the test particle. When the test particle is far from the planet it 
follows a roughly Keplerian orbit around the star. When the test particle is very close to the 
planet, it also follows a roughly Keplerian orbit, but this time around the planet. While it is 
close to the planet, the test particle experiences a tidal force from the star, with a magnitude 
given by 

F. ^ 2MpL^, (49) 
a^ ap 

where r is the distance between the test particle and the planet, G is the gravitational 
constant, is the mass of the central star, etp is the semimajor axis of the planet's orbit, 
and nig is the mass of the dust grain. The Hill radius is the distance r at which this tidal 
force equals the gravitational force of the planet on the test particle, GMprUg/r^, or 

^-n^J a„ (50) 

where the factor of 3 is included for historical reasons. 

We make use of the Opik approximation; while the test body is inside the Hill sphere of 
the planet, we assume that the motion is described by the two body problem in the frame 
revolving with the planet. This involves ignoring both the tidal force of the star and the 
Coriolis force associated with the motion of the planet around the star. The acceleration 
due to the Coriohs force is 

where fip is the mean motion of the planet (the angular speed of the planet's revolution 
about the star) and V is the velocity of the test particle in the rotating frame. A rough 
estimate for V is v^,^ the Keplerian velocity of the planet, so the Coriolis force is larger than 

— 1/3 

the tidal force by a factor Oyp/rn ^ l^p at the Hill radius; in this section jjip = Mp/M^. 



4.1.1. Inside the Hill sphere 

We assume that the test particle orbits the star rather than the planet. Hence the 
particle follows a hyperbolic orbit relative to the planet during the close encounter. The 
hyperbola is specified by its eccentricity e and the transverse semimajor axis a (note that 
the semimajor axis of the dust particle a should not be confused with the radius a of the 
dust particle). Barred elements arc calculated relative to the planet. Note that e > 1 and 
a < 0. Alternately we can specify the specific energy E — —GMp/2'h and specific angular 
momentum L — -yCMj, a(l — e^) relative to the planet. The energy is positive, so the 
particle would have a finite velocity as it traveled to infinity, if we ignore the effects of the 
star. This velocity is known as the velocity at infinity; we denote it by U . 
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The approximations described above allow us to find an analytic relation between the 
angle of deflection resulting from the encounter and the eccentricity of the (hyperbolic) orbit 
of the test body around the planet. We then relate the eccentricity to the periapse distance 
q between the planet and the test particle, and the velocity at infinity. 

The specific energy and angular momentum of a test mass interacting with a planet of 
mass Mp are given by 

^=V-^ (52) 

2 f ^ ' 

and 

,63) 

In these expressions G is the gravitational constant, V is the velocity of the test mass relative 
to the planet, and f is the distance between the two bodies. The angle 9 is measured from 
a fixed line (which we choose to be the apsidal line, which connects the star and the planet 
at the point of closest approach) and the line joining the two bodies. 

The distance f is given by 

f=- ^ -. (54) 

1 + e cos 61 ^ ^ 

The quantity p = L'^/GMp is called the semilatus rectum. The periapse distance (at 
closest approach) is denoted by q and is given by g = p/(l + e). 

The semilatus rectum p is related to e and a by 

p = a(l - e^) (55) 

Note that p > 0, and that q = a(l — e)>0. Using the definition of a, we can write the 
specific energy in terms of p and e: 

E=^(e-'~l). (56) 
2p 



Prom eqn. (52) 



E=-v^,-^ = -U\ (67) 



where Vg is the speed at closest approach, and U is the speed at infinity. 
Prom equation (54), as r — > oo, the angle 9 tends to 



^0 = — arccos ^ (58) 
e 
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Referring to Figure (5), and denoting the angle of deflection by 7, simple geometry gives 

7 TT 1 

sin - = sin(6lo - -) = - cos 6^0 = |. (59) 
Since e = 1 — {q/ a), we find 

^ ^ ■ ^ ^ (60) 



sm 



1 + 



GMr, 



It is convenient to write this in terms of the impact parameter s. From conservation of 
angular momentum, sU = qVq. We find 



172" 



Using this in eqn. (60) we find 



. 7 
sm — 
2 



4„2 



1 + 



GMp 

w 



-1/2 



(61) 



(62) 



4.1.2. Relating U to the orbital elements a, and e of the test particle 

We will assume that the impact parameter s is uniformly distributed between and 
Tij. However, we still need to know U, the speed at infinity, in the frame rotating with the 
planet. We can find U in terms of the semimajor axis a and eccentricity e of the test particle, 
where the orbital elements are now calculated relative to the star. The relations between E, 

L, a, e, and p are as given above, but with the star playing the role of the central mass, so 
that Mp is replaced by M^, the mass of the star. Since the test particle is bound to the star 
E < 0, a > 0, and e < 1, at least until the scattering event that ejects the test particle. 
We will restrict our attention to the case where the planet and the test particle orbit in the 
same plane. 

We start by dividing the particle velocity (relative to an inertial frame centered on the 
center of mass) into a radial Vr and a transverse Vt part. From the definition of angular 
momentum, the transverse velocity at the time of the close encounter is 

Vt = L/r ^ vpj—{l-e^), (63) 

where Vp — yjGM^j ap is the Keplerian velocity of the planet, and we have made use of the 
fact that the star-particle distance r ^p (it can differ by rjj) during the encounter. 



-23- 



The magnitude of the total velocity v is found from the expression for the energy, 



E = -V 

2 Qjp 



2a ' 



(64) 



where we have again set r = in the first equality. Solving equations (63) and (64) for Vr, 

1/2 



2-^ 

a 



(1 - e^) 



(65) 



Note that if the orbits just cross, a(l — e)/ap = l, l + e = 2— a^/ a, and Vj. = 0; the collision 
occurs at periapse. Then U'^ = [vt — ~ v^. However, the orbit of the test particle need 
only pass through the Hill sphere, so that a(l — e)/ap ^ 1 ± r///ap ^ 1 ± /ip^'^, and the 
radial velocity Vr ~ ^/2fip^'vp, which for a Jupiter mass planet is about 0A5vp. If there are 
multiple planets in the system, the test particle periapse need not be comparable to ap, in 
which case Vr could be slightly larger than Vp. The transverse velocity Vf is always of order 



Up. 



To find U, we transform to the planet frame, which entails subtracting Vp from vf. 

Ur^Vr ^ V^/^yS (66) 

(67) 



Ut^vt 



Vp ~ 



dp 



2nl'^ + (3 - 2V2) 



(68) 



4.. 1.3. The change in energy 

In the frame rotating with the planet, the close encounter simply rotates the test particle 
velocity by the angle 7, without changing its magnitude U . The components of the scattered 
velocity are 

— f/r cos 7 — f/t sin 7 (69) 
U[ = [/^sin7 + t/t C0S7. (70) 

Transforming back to the inertial frame, the square of the new velocity is given by 

^/2 ^ ^/2 ^ ^ ^ jj2 ^ jj2 ^ 2Vp{Ur siu J + Ut COS 7) + V^. (71) 

The change in specific energy 

AE = ^ ^ « Vp [Ur sin 7 + Ut{l - cos 7)] (72) 
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The radius r does not change appreciably during the encounter. Using the relations between 
U and V, 

AE — Vp [vr sin ^ + {vt — Vp){l — cos 7)] . (73) 

4.1.4- The small scattering angle approximation 

Equation (68) shows that C/^ O.AVp for test particles with a >> a^ scattering off 
Jupiter mass planets. Prom equation (62), we find 

7 « 2GMp/{sU^) = 2^/.f (^) . (74) 

The small angle approximation is valid for s > s^, where 

S7 ~ /^p (^'^^ ^'^^^ 

There is another constraint on the minimum value of the impact parameter s, namely 
that the test particle does not physically collide with the planet. This defines Smin, the 
impact parameter for wliidi tlie i)eiiapse q — rp, where is the radius of the planet. For 
a >> ap, Smin ~ rp-\/2/Xp(wp/?7)2 s-p/fp. At 5AU s^j„ ~ 7rp for a planet with Jupiter's mass 
and radius. Thus non-coUisional close encounters occur if the impact parameter satisfies 

7rp < s < th, (76) 

For a Jupiter mass planet at 5 AU, s-y ~ 25rp, about four times larger than the minimum 
impact parameter. For the same values of Mp and ap, th ~ 750rp. In other words, the 
small angle approximation is valid for almost any impact parameter that does not lead to a 
physical collision. 

Assuming 7 << 1, the change in energy due to a single close encounter for a test particle 
with a >> ap is 

A£;(/.p, ap, .) ^ 2'l-r^i,f (^)' (^) vl (77) 

4.1.5. The ejection velocity 

We are now in position to estimate the typical ejection velocity of a test particle interact- 
ing with a Jupiter mass planet. The particle must have an energy satisfying —A.E < E < 
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when the last encounter occurs. On average it will emerge with an energy E AE/2 after 
the encounter, so it is ejected with a velocity 

vooil^,, a„ s) ^ 24/'' (^) ^v,. (78) 

Note that the ejection velocity scales as the 1/2 power of the planet's semimajor axis Bp. 
The scaling with the mass ratio /ip is complicated by the appearance of the Safronov number 
U/vp. Prom equation (68), the Safronov number is independent of /ip for /Xj, << 6 x 10~^, 
but scales for larger [Xp. Hence the ejection velocity will scale for Saturn 

mass planets, and as ji]/'^ for planets significantly more massive than Jupiter. If the planet 
has a substantial eccentricity, as many extrasolar planets do, then Vt ~ Vp and the ejection 
velocity will scale 

For a Jupiter mass planet at 5 AU from a solar mass star, 

Let X — s/th] then the cross section for an encounter with impact parameter between s 
and s + is 2t: sds — 27rrfj xdx. We average over x to find the mean escape velocity. We 
integrate from x^ = s^/th ~ fJ^p^^ << 1 to s = Oirn, where a is a dimensionless constant of 
order unity, and is the minimum value of s for which equation (74) is valid. Doing the 
integration, 

I nxdx 6 J a 

where we have neglected terms of order fJ^pVr^{rH)', we use the notation Voo{rH) — v^oil^p, ^p, s — 
th)- The rms escape velocity is calculated in a similar manner. We find 

-Voo^Th), (81) 

a 

so the spread in escape velocities is similar to the mean escape velocity. 

For those rare encounters with Smin < s < s^, the scattering angle is of order unity, and 
the velocity kick experienced by the particle is larger than the estimate in eqn. (78). We 
take sin 7 f» cos 7 ~ 1/ -\/2, so 

AE^vl, (82) 

and 

^"00 ~ Vp, (83) 

much larger than the small angle limit. About 1% of the particles will be ejected with 
velocities of tens of kilometers per second. 
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4.1.6. Ejection time scale compared to collisional timescale 



So far wc have assumed that the grains are ejected before they suffer sufficient colUsions 
with other dust particles to alter their orbits substantially. We are now in a position to 
check this assumption. We start with an estimate of the dust-dust collision time. The 
optical depth r for the most massive debris disks is of order 0.001. This optical depth is due 
to particles with radii a ~ A/27r, where A is the wavelength at which the disk is observed. 
Since the disks are usually detected by their IR excesses, the dust particles responsible for 
the optical depth have a ~ 0.3//m or smaller. We are interested in larger ejected particles, 
with fleject = 25/xm. These larger particles will suffer a collisions with a particle of size 
o-target after roughly 1/T{atarget) orbits, where T{atarget) is the optical depth to particles of 

size a^arget- 

Using the scaling mfm oc m~^, the optical depth in particles of size atarget, as viewed at 
a wavelength of order or smaller than atarget, is 



The typical number of orbits required for our test (ejected) particles to sweep up their 
own mass in smaller particles, and hence have their orbits altered substantially, or to strike 
a larger particle, is 



orbital periods. 

Thus the most efficient way to alter the orbit of a particle is to collide with another 
particle of the same size, assuming that the mass is logarithmically distributed, and assuming 
that gravity plays no role in the collision. The typical number of orbits for a 25//m size 
particle to coUide with a particle of its own mass is ~ ^/^{a — 25//m) 7 x 10^. 

Dust-dust collisions can also shatter the grains. The binding energy for a target grain 
of mass m is (m. A collision with a grain of mass nic moving with relative velocity v will 
shatter the target if the kinetic energy exceeds the binding energy of the target. Allowing 
for the possibility that only a fraction fxE of the kinetic energy is available to disrupt the 
target, the critical mass needed to shatter the target is 




(84) 




(85) 



(86) 
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In deriving this result we assume that the relative velocity is evp, corresponding to the 
random velocity of a grain suffering close encounters with a massive planet of semimajor 
clXlS dip . The eccentricity of the target grain undergoes a random walk in e, starting at zero 
and evolving to e = 1, e{t) ~ \/Det-, so wc take e = 2/3. The minimum mass needed to 
shatter the target grain is then rric ~ 0.5( ap/5AU)m//xE, roughly a mass equal to that of 
the target. 

How long does it take for the typical particle to be ejected? The test particles undergo 
a random walk in energy, with a step size given by equation (77). The number of collisions 
needed to random walk from E = -GM /2&p = to E = is iV a; S^/^/lG/Xp^l The 

probability of a close encounter on each periapse passage is 

(87) 



{2ti h-pTh) 2 V 3 
The total number of orbits up to ejection is 

3 . /10-3N^/^ 



N.,ect ^ « 4 X 10^ f HZ^ 



We conclude that particles larger than about 25/im will be ejected before their orbits are 
significantly altered by collisions with smaller or similar size dust grains. 



4.1.7. Numerical Results 

We have carried out numerical integrations of test particles in the gravitational field of 
one or more massive planets orbiting a star. We use the publicly available SWIFT (Levison 
& Duncan 1994) integration package, which is based on the Wisdom & Holman (1991) sym- 
plectic integration scheme. Wc used the rmvs3 integrator, in order to integrate through close 
encounters between the test particles and the planets. All integrations were for 100 million 
years, with a timestep chosen to be small enough to resolve periapse passage for the esti- 
mated most extreme test particle orbits, or to resolve the perijove passage for planet grazing 
close encounters (taken to be at 2rp) , whichever was smaller. If a particle passed within 2rp 
is was deemed to have collided with the planet and was removed from the integration. 

The initial eccentricities and inclinations (measured from the planet's orbital plane) of 
the test particles were generally set to 0.1 and 0.087 radians respectively, although we tried 
runs with other values. The final ejection velocities did not depend strongly on the initial 
e and i. The test particles were given semimajor axes ranging from 0.5 to 1.5 a^ in single 
planet cases, since in those cases particles starting at larger distances from the planet were 
typically not ejected by the time the integrations were halted. 
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Figure (6) shows the result of a numerical integration of ~ 600 test particles in the 
gravitational field of Jupiter and the Sun. Our analytic calculations should be a good 
approximation to this case, since ej = 0.048. The test particles were started with a range of 
semimajor axes between x and y, with initial eccentricities e ~ 0.05 and inclinations i ~ 0.05 
radians. The figure shows the velocities at which the particles were ejected, including the 
(small) correction for their finite distance from the Sun when the integration was stopped. 
The stopping criterion was that the test particle have a positive energy relative to the Sun 
and a semimajor axis larger than lOOAU. The mean escape velocity is ~ Ikms"-*^, shghtly 
lower than predicted by equation (79). 

We found that the final inclinations were small, typically within 10 degrees of the orbital 
plane of the planet, but with a significant minority of particles ejected at inclinations up to 
30 degrees. 

We checked the scaling of Vej with a^ by varying the size of the planet's orbit over a 
decade. A least squares fit to the mean ejection velocity as a function of a^ gives Vgj ~ 
with an error of about 0.05 in the exponent. This is consistent with the predicted slope of 
—0.5 given by equation (78). We also checked the scahng with Mp, finding Vej ~ M°-^^, 
compared with the prediction exponent of 1/4. In summary, equations (78) and (79) appear 
to be a good description of the ejection process. 

The numerical integrations employing multiple Jupiter mass planets gave similar results, 
with the test particle ejection velocity determined by the most massive of the planets in the 
integration. The same was true of integrations involving a single Jupiter mass planet together 
with several 1-10 Earth mass bodies. 

These results suggest that equations (79) and (81) can be used to describe the ejection 
velocity of small particles in systems with one or more Jupiter mass planets. 

4.2. Dust Luminosity 

In order to estimate the (age-dependent) dust luminosity of a Vega-like star, we must 
first estimate the number of grains in the disk. Spangler et al. (2001) have measured IR 
excesses for stars in several nearby clusters. Their results are consistent with the following 
relation: dust disk mass oc (age)~^. We normahze this relation by estimating the (size- 
dependent) number of grains in the disk around /3 Pictoris, with an age ~ 12 Myr (Zuckerman 
et al. 2001). 

We assume that the grains around P Pic are located at a distance D — lOOAU from 
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the star and that they have a size distribution dN/da — Ca~^-^, where N{a) is the number 
of grains with size < a. This is the equihbrium size distribution that results when mass is 
redistributed among grain sizes via shattering in colhsions (Dohnanyi 1969). We obtain the 
normahzation constant C by setting the emission at 800/im equal to the observed value of 
Fy = (115 ± 30)mJy (Zuckerman & Bccklin 1993). The specific fiux is given by 



where the absorption cross section is Qabs^ra^, Bj^{T) is the Planck function, and the grain 
temperature T(a) is determined by equating the absorption and emission rates. For /3 Pic, 
we take effective temperature T^s — 8200 K, luminosity L — S.7Lq, and radius R — 1A7Rq 
(Crifo et al. 1997). Grains with a > 1cm contribute little to F,^(800/xm). Terminating the 
integral in equation (89) at a = 1 cm yields C = 4.65 x 10^^ cm^-^ for sihcate grains, implying 
a disk mass of 1.4 x 10^^ g (~ 19 lunar masses, 0.23 Earth masses) in grains with a < 1cm. 
If we assume graphite (water ice) rather than silicate composition, then the estimated mass 
increases by a factor of 1.6 (4.6). Since pure water ice grains are very unlikely, our mass 
estimate is robust against variations in compositions. If the grains are located 30 pc (500 pc) 
rather than 100 pc from the star, then the estimated mass decreases by a factor of 2 (increases 
by a factor of 3). Thus, though the dust actually will be distributed over a range of distances 
from the star, this does not seriously affect the estimate of the total dust mass. 

Our result can be compared with the analysis of Li & Grccnbcrg (1998), who constructed 
a detailed dust model to account for the emission (at multiple wavelengths) from the dust 
around P Pic, using fiuffy grains. In their model, the total dust mass needed to produce the 
observed emission is ~ 2 x 10^^ g, lying mostly beyond 100 AU from the star. 

It is interesting to note that estimates of the dust mass in circumstellar disks typically 
assume that the grains are in the Rayleigh limit, i.e., that a -C A, the wavelength of the 
observation (see, e.g., the review by Zuckerman 2001). In this regime, Qabs oc a. In the 
opposite regime (a ^ A), Qabs ~ 1- Suppose oq is the minimum a for which Qabs ~ 1, for 
fixed A. Then, Qabs ~ ct/oo in the Rayleigh limit. The opacity k (absorption cross section 
per unit mass) of the emitting grains is the crucial quantity for estimating dust mass from 
emission; k = 3(5abs/4ap. In the Rayleigh limit, k ~ 3/4aop. In our estimate of the dust mass 
(with our assumption that dN/da oc a~^'^), the emission is dominated by grains with a ~ aoi 
since these arc the most abundant grains with Qabs ~ 1- Thus, k ~ 3/4aoP for our estimate 
as well as for the Rayleigh limit, and these two different methods yield approximately the 
same mass in grains contributing significantly to the observed emission. We feel that our 
scenario, with a Dohnanyi size distribution, is more realistic than the standard scenario, in 
which all of the grains that contribute significantly to the emission are supposed to be in the 
Rayleigh hmit. 




(89) 
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Combining the above estimate for the number of grains in the disk around /3 Pic with 
the Spangler et al. (2001) result, we estimate 



, ~ 4.7 X 10^^ cm^-^ — a--^-^ , a < 1 cm , (90) 

da \ t / 

where t is the age of the system. 

The Spangler et al. (2001) result suggests that grains are lost at a rate dN/dt ~ —2N/t. 
The following could be important sinks for the grains: 1. gravitational ejection from the 
system, 2. gravitational ejection into the star, 3. incorporation into planets or planetesimals, 
and 4. ejection of small shattering fragments by radiation pressure. We will not attempt to 
estimate the relative importance of these sinks in this paper. Rather, we will introduce the 
unknown factor /ej, equal to the fraction of the grains that are lost by gravitational ejection. 
Then, the specific dust luminosity is 

L.At,Vei,a) = 7.8 X 10^«yr-icm-i(cm/s)-i (1^]^^ /^^/^(^^.) , (91) 

where we have taken time equal to the age of the star when the grains are ejected and 
fy{vei)dvej is the fraction of the grains ejected with speed between Vej and Vej + dvey"^ 

At very early ages, the primordial gas disk will not yet have dissipated, and the drag 
force on grains could dramatically reduce the dust luminosity. Also, we do not expect large 
grains to be ejected prior to planet formation. Thus, L^^a = for i < tcr, an unknown critical 
age. We will assume that tcr lies in the range 3 to 10 Myr. 



4.3. Flux at Earth 

Assuming /survive — 1) the flux, at Earth, of grains with size > a, F{a), is given by 

F(a) = 8.2xlO-5yr-ikm-Vbeam/ej (j^) / d^ei fv{v,i)'^Q{U ^ t-d/v,^) , 

v^, ^^^^ 



where 

) ^ej ^cr 



^We assume that fy is independent of a. 
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If we assume Vd^e/vej ~ 10, this yields a flux of 2.6xl0~^^(25//m/a)^-^(10pc/d)^cm~^ s~^, 
which may be compared with Fig. (1); the AMOR point source supphes a flux at Earth of 
10~^^cm~^ s~^. It would appear that the point source is not due to a debris disk. 

We will be interested in sources within the LB, so /survive = 1 when d < tb- With 
$ ^ 0.94 V and p — 3.5gcm~^ (appropriate for sihcate grains in the LB), equation (45) 
yields a critical grain radius 

/ ^ N 1/2 / \ -1/2 

/survive = 1 whcH a > a^r- Wc will usc v^^^q as a good approximation to f^,© in equation (94). 

Suppose the threshold of detectability requires 20 events per year, so that the threshold 
flux Fth = 20yr-M-|. Then, we can estimate the threshold distance dth out to which 
sources might be detected by making the optimistic assumptions that tcj = ^cr and /ej = 1 
(we will also assume that /beam = !)• Inserting equation (94) into equation (92) and assuming 
Vd,e/vei ~ 10 yields 

4 \ '1/13 / , \ -12/13 

-^COl » I f^r 



In Table 3 we give d^h and a^r for the cases that Tcr = 3 or lOMyr and Acoi — 10^ km^ or 
10^ km^; recaU that AMOR has A^i fa 10 km^ 

Since the flux drops off as i^^, it is highly advantageous to catch the star when tej ~ ta. 
In Figure 4, we plot, as a function of d, the upper and lower ages for which tcr < ^ej ^ '^tcr, 
assuming a range of ejection speeds between 0.5kms~^ and 2kms~^. Since the total volume 
and favorable age range both increase with d, distant stars that are possibly too old to have 
observable IR excesses today may be a significant source population, despite the d^"^ decrease 
in flux. However, nearby identifiable Vega-like stars are more attractive sources, since we 
could then combine the dust flux information with other observations to learn more about 
a particular object. 



4.4. Candidate Sources 

4.4- 1- Gliese Catalog Stars 

Are there any nearby Vega-like stars for which we can expect a detectable dust flux at 
Earth? To answer this question, we consider stars from the Ghese Catalog with far-infrared 
excesses from IRAS (see Table X in Backman & Paresce 1993). 
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In Table 4, we give each star's distance d and its space velocity {U, V, W) with respect 
to the Sun. In calculating these quantities, we have taken coordinates, parallaxes, proper 
motions, and radial velocities from SIMBAD.* We estimate a star's age T from its IR excess 
r, which is the fraction of the star's luminosity that is re-emitted in the IR by dust. Adopting 
the Spangler et al. (2001) result that r oc t^'^, we take 



t = t(/3Pic) 



>(/3Pic) 



0.5 



(96) 



The dust fluxes in Table 4 were calculated using equation (92) with limiting grain 
size a = max(10/im, a^) (ocr is taken from eq. 94 with v = v*,©), a flat distribution of 
ejection velocities between 0.5 km s""*^ and Skms^"*^, tcr = 3Myr, and /beam = /cj = 1- Thus, 
these estimates are optimistic. In §7 below, we suggest how to build a radar detector with 
Acoi ~ 3 X 10^ km^. Four of the stars in Table 4 would yield 20 or more meteors per year 
(i.e., F > 6 X 10~^yr~^ km~^) for such a system. Since igj is substantially greater than 
for most of the stars, increasing tcr to 10 Myr yields reduced fluxes for only three stars: Gl 
219, 297.1, and 673.1. The flux for Gl 219 (/3Pic) vanishes, whereas the fluxes for Gl 297.1 
and 673.1 remain greater than 2 x 10~^ yr~^ km~^. 

The apparent location of a dust source on the sky is determined by the velocity of the 
dust relative to the Sun (eq. 4). If the dust ejection speed Vej were much greater than the 
star's speed v^^q, then the dust would appear to come from the location of the star itself 
(since the speed of light is much greater than the star's speed). However, v^^ is typically an 
order of magnitude or more smaller than w*,©. Thus, the location of the dust source on the 
sky is primarily determined by v*,© and need not be anywhere near the actual location of 
the source. 

Since -u^, © is determined in part by the velocity %,,lsr of the Sun relative to the local 
standard of rest (LSR), we may expect some degree of clustering of the apparent directions 
around the solar apex (i.e., the direction of '%,lsr)- Dchncn & Binncy (1998) give U = 
10.00 ± 0.36, V = 5.25 ± 0.62, and W = 7.17 ± 0.38 kms"^ for u©,lsr- They find that 
the velocity of the Sun with respect to nearby young stars is similar, except that V ~ 
12kms~^ in this case. The coordinates of the solar apex with respect to the LSR are 
(Z,6) = (27.7°, 32.4°) or (A,/3) = (248.4°, 32.2°), while those with respect to nearby young 
stars are (/, b) = (50.2°, 24.7°) or (A, /3) = (265.7°, 48.9°). In addition to Galactic coordinates 
{l,b), we also make use of ecliptic coordinates: ecliptic longitude A and ecliptic latitude (3. 



^Note that we exclude one star from Table 4 in Backman & Paresce (1993), Gl 245, because SIMBAD 
does not give its radial velocity. 
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These coordinates are useful in radar studies of meteors, because most of the observed grains 
originate in the echptic plane. In Table 4, we give the apparent direction to the dust stream 
in ecliptic coordinates. Figure 7 is a plot of these directions on the sky. 

4-4- 2- Young Clusters 

Nearby young clusters could potentially yield strong dust fluxes. The cluster members 
may not have appeared in searches for the Vega Phenomenon, since the stellar luminosities 

may be too low. However, most nearby clusters have poor tuning between distance and age. 
There are some very young clusters that are too far away for the dust to have reached us 
yet and some nearer clusters that are much older (see Table 1 in Spangler et al. 2001). 

The recently discovered Tucana Association (Zuckerman, Song, & Webb 2001) might be 
more suitable. Zuckerman & Webb (2000) find a distance ~ 45 pc, age ~ 40 Myr and space 
velocity ([/, V, W) ^ (-11, -21, 0) kms"^ Stelzer & Neuhauser (2000) find a younger age of 
10-30 Myr. Adopting numbers from Zuckerman & Webb and making use of the convenient 
coincidence that 1 kms~^ = 1.02 pc/ Myr, we find that dust ejected at speed fcj = 1.5kms~"^ 
when the stars were 15 Myr old would be reaching us today. The position of the dust stream 
on the sky would be (A, f3) = (126°, -46°). 

The Pleiades cluster is another potential source, with a distance of 118 pc and an age 
of ~ 120 Myr (Spangler et al. 2001). Although the large distance would suppress the flux 
(because of both the dependence in eq. 92 and the need for larger grains in order to 
avoid magnetic deflection) , this could be partially compensated by the large number of stars 
in the cluster. Robichon et al (1999) flnd {U,V,W) = (-6.35, -24.37, -13.02) km s~^ for 
the cluster motion. Ignoring the Galactic potential, this would yield (A, /3) (276°, 73°) for 
the position of the dust stream. 

4.5. Has Dust from /3Pic Been Detected? 

Baggaley (2000) detected a "discrete" source of radar meteoroids with a central location 
(A, /3) (280°, —56°), an angular diameter of ~ 30°, and a dust speed relative to the Sun 
Vd,Q ~ 13kms~^. He claimed that the central location and dust speed could be reproduced 
if pPic were the source and the grains were ejected with speed Vej ~ 29kms~^. 

This result seems unlikely on theoretical grounds, since it is not clear how the grains can 
be ejected with such large speeds. In addition, we do not flnd that dust emitted from /3 Pic 
reproduces the location of the discrete source for any v^j ■ Baggaley did his calculations in the 
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LSR frame, taking the Sun's motion with respect to the LSR from Binney & Tremaine (1987), 

who give {Uq, Vq, Wq) = (9, 12, 7) km s~^. Baggaley found the direction to the discrete source 
to be = (49°, -72°), (Lb) = (259°, -28°), whereas we find (A,/5) = (58.6°, -81.3°), 

(/, 6) = (267.4°, —34.1°). Note that the two coordinate pairs given by Baggaley do not 
actually correspond to the same point on the sky 

In Figure 8, we plot the sky position of a dust stream emitted by /?Pic for various values 
of Vcj (triangles). The box indicates the central location of Baggaley's discrete source. For 
reasonable ejection speeds (fcj < 3kms^^), the location of the dust stream differs dramati- 
cally from that of the discrete source. For Ucj ~ 30kms~^, the locations arc much closer and 
Vd,Q ~ 13kms~\ Although the dust stream does lie within the ~ 30° extent of the source 
in this case, it is still 12° away from the center of the source. 

The particle flux coming from the discrete source is several orders of magnitude larger 
than we would expect from a debris disk at 20 parsecs. From eqn. (92) we find that Beta 
Pic produces a flux of a = 10/xm particles of 10~^ per square kilometer per year at Earth. 
Given our estimate of the collecting area of the AMOR detector, somewhat less than 10 km^, 
AMOR is not capable of detecting such a low flux. 

Thus, it appears that Baggaley's discrete source is not related to (3 Pic. 

4.6. Potential Distributed Sources 

In addition to the discrete source. Figures 2b and c in Baggaley (2000) suggest the 
presence of a distributed, band-like feature. Here we discuss three possible sources for a 
distributed feature, namely the Galactic plane, Gould's Belt, and the spiraling of grains in 
the local magnetic field. 

4-6.1. The Galactic Plane 

The scale height of young stars above the Galactic plane is ~ 90 pc (Gilmore & Reid 
1983). Thus, we might expect to see a signature of the plane in the dust fiux since we can 
detect grains from sources beyond 100 pc if a and v are large enough (but v must not be 
so large that the grains are destroyed). The effects of the Galactic gravitational potential 
cannot be ignored for such large distances, but we do not expect the Galactic potential to 
deflect grains out of the plane. As mentioned in §4.4.1, the apparent direction of a dust 
stream with a low ejection speed is primarily determined by the velocity of the source with 
respect to the Sun. Since the Galactic plane is a feature in physical space rather than velocity 
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space, it may seem that it should not actually appear as a band on the sky in a dust flux 
map. However, if the radial velocity of a source (with respect to the Sun) exceeds the dust 
ejection speed, then, unless the source passed near the current location of the Sun at some 
time in the past, the dust will never reach the Sun. Thus, though the velocity of the source 
with respect to the Sun determines the direction of the dust stream, the location of the 
source does restrict which velocities can give rise to an observable dust stream. 

From equations (92) and (94), the dust flux F oc c?"^'^^, so it is not clear that distant 
stars in the plane can indeed produce a noticeable band on top of the relatively isotropic 
distribution produced by nearby stars. The following simple estimate shows that a band is 
expected, although at a low contrast level. For stars with distance i? < 90 pc, the volume 
element (in which young stars are found) is dV — AnR'^dR, while dV — {90pc)2TTRdR for 
i? > 90 pc. Thus, the ratio of the flux Fpiane due to the Galactic plane to the isotropic flux 
Fiso (due to nearby stars) is 



(180pc)7r C dRR-^-^^ 

plane \ ' J 90 pc 



= 1.6 , (97) 



where we have multiplied each integrand by a factor R to account for the increased probability 
of catching a star when Tej T^. 

Since the components of random stellar velocities perpendicular to the plane are smaller 
than those parallel to the plane, the Galactic plane should appear as a fuzzy band on the 
sky even considering only stars closer than 90 pc. 

Note that the observed dust flux band should be warped in such a way as to favor the 
direction of the solar apex. 



4.6.2. Gould's Belt 

Gould's Belt is a band on the sky, inclined ^ 20° with respect to the Galactic plane 
and with an ascending node ~ 280°, along which young stars tend to lie (see Poppel 1997 
and de Zeeuw et al. 2001 for reviews). From a Hipparcos study of early-type stars, Torra, 
Fernandez, & Figueras (2000) concluded that fa 60% of the stars younger than 60Myr 
and within 600 pc of the Sun are in Gould's Belt. OB associations, young star clusters, and 
molecular clouds have long been known to trace Gould's Belt. Recently, Guillout et al. (1998) 
detected late-type stellar members of the Gould Belt population, by cross-correlating the 
ROSAT All-Sky Survey with the Tycho catalog. Thus, Gould's Belt could be responsible for 
a distributed grain flux feature. 
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Guillout et al. (1998) found that the three-dimensional structure of Gould's Belt is disk- 
like. The outer rim of the Gould Disk is ellipsoidal, with a semi-major axis of ~ 500 pc and 
a semi-minor axis of ~ 340 pc. The center of the structure is ~ 200 pc from the Sun, towards 
/ ^ 130°. On the near side (/ ^ 310°), the inner edge of the Gould Disk lies only fti 30 pc 
from the Sun, although most of the young stars lie beyond ~ 80 pc. 

Estimates of the age of Gould's Belt range from 20-90 Myr (see Torra et al. 2000 
for a review); Torra et al. favor an age between 30 and 60 Myr. Grains traveling with 
fcj = 1.5kms~^ reach a distance of ^ 90 pc from their star of origin in 60 Myr. Thus, only 
the near side of the Gould Disk should be visible in the dust flux. 

The kinematics of Gould's Belt is still poorly understood, but it is clear that the mem- 
ber stars are undergoing some sort of expansion. If the expansion dominates the random 
velocities, then the near side of Gould's Belt should appear on the opposite side of the sky 
in the dust flux, i.e., centered on / pa 130°. If the random velocities are more important, 
then we would expect to see a feature lying in the same direction as the near side of Gould's 
Belt, due to stars moving towards us. 

4-6.3. Spiraling of Grains in the Local Magnetic Field 

Thus far, we have only considered grains that are not dramatically deflected as they 
travel from their source to the Earth. Suppose the magnetic fleld were uniform throughout 
the Local Bubble. In this case, deflected grains would produce a wide band centered on the 
plane perpendicular to the field. The actual field presumably includes a significant random 
component, so the resulting distributed feature (if any) may be more complicated than a 
wide band. 



5. DUST FROM AGB STARS 

Stars with main sequence mass M < 6Mq spend time on the asymptotic giant branch 
(AGB). During this evolutionary phase, a wind is driven off the star, with mass loss rates as 
high as ~ lO~^M0yr~^. Grain formation in these outfiows hkely supphes a large fraction of 
the dust in the ISM. See Habing (1996) for a general review of AGB stars. 

Since the grains spend only a limited time in a high-density environment, the size dis- 
tribution of the newly formed dust is thought to be dominated by grains with a < l/im (see, 
e.g., Kriiger & Sedlmayr 1997). Observations confirm the significant presence of submicron 
grains. Jura (1996) examined the extinction due to circumstellar dust around seven 0-rich 
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giants and found that it rises towards the ultraviolet. Prom observations of scattering, Groe- 
newegen (1997) found a mean grain size a ~ 0.16//m for the circumstellar envelope of the 
carbon star IRC +10 216. 

Even if submicron grains dominate the outflows from AGB stars, the luminosity of 
grains with a > 10//m may be high enough to result in detectable fluxes at Earth, since the 
dust production rate is so high. Grains with isotopic ratios indicative of formation in AGB 
stars have been discovered in meteorites. Although most of these grains have a ~ 0.5/im, 
some are as large as a a; 10/xm (Zinner 1998). 

Suppose the grains all have radius a and that a total dust mass Md is emitted during the 
AGB lifetime r. The dust luminosity is then L = ?,Md/{Anpa''T). We adopt ~ O.OIMq 
and the canonical lifetime r ~ 10^ yr (Habing 1996). The flux at Earth follows from equation 
(3). For simplicity, we take /beam = 1 and /survive = 1; also, since the typical grain ejection 
speed in an AGB wind is ~ lOkms^^ (Habing 1996), we take Vd,Q/vej = 1. Thus, the flux 
as a function of AGB star distance d is approximately 

F{d) ~ 114 yr km {^^^-2^^ (3.5 gcm-^) (lO^im) (lO^ yr) (lOOpc) 

(98) 

In order for the grains to be traceable to their point of origin, we require that the 
gyroradius tb > d. This yields a minimum acceptable grain size (eq. 45), so that 

^' ^ ^ V 10-^^0 / V 3-5 gem- V VO-SVy \^liG ) 

Viokms-V vio^yv Vioopcy ^ ' 

The above flux is for the optimistic case that all of the dust mass is in the smallest grain size 
at (25//m in this case) such that the grain will barely be undeflected. In cgs units, this flux 
is 2.5 X 10~^^cm-^s~^ tantahzingly close to the flux of the point source seen by Baggaley 
(2001). 

Jackson, Ivezic, & Knapp (2002) have studied the distribution of AGB stars with mass 
loss rates in the range lO^^-lO^^M© yr"^. They find that the number density of AGB stars 
TiAGB ~ 4.4 X 10~^ pc~^ in the solar neighborhood. Thus, the number of sources within a 
distance d from the Sun is N{d) ^ 1.84(d/100pc)^. 

In Table 5 we give the flux F from equation (99) and the number N of sources for a 
few values of d. We also give the threshold grain size a* and the fraction g of the emitted 
dust mass that would need to be in such large grains in order to yield an observed flux at 
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Earth of 2 x 10 ^yr ^km ^ (corresponding to 20 events per year with a collecting area of 
10^ W). We see from the table that AGB stars could be a significant source population for 
radar detection of extra-solar meteoroids if a small fraction of the grains in the outflow are 
large. Note that the source AGB star will not be visible as such, due to the large distances 
to AGB stars and their short lifetimes. Also, the distances are large enough that the grain 
trajectories will be affected by the Galactic potential. 

Even if all the grains in the AGB outflow are smaller than lO^m, there is still a chance 
to see the AGB star as a source on the sky, since some AGB stars apparently have long-lived 
disks in which grain coagulation occurs. This phenomenon has been observed for binary 
stars; the presence of the companion apparently causes some of the outflowing grains to be 
deflected into a disk. A well-studied example is AC Her (Jura, Chen, & Werner 2000). Jura 
et al. find that the grains must have a > 200/xm in order to remain gravitationally bound, 
due to the extreme luminosity of the evolved star. They estimate the mass of the dust disk 
(in grains with a < 0.1 cm) to be ~ 1.2 x 10^° g. Their observations of infrared emission imply 
the presence of much smaller grains as well. They suggest that shattering during collisions 
between disk grains leads to a radiation pressure-driven outflow of smaller grains. The 
lifetime of such disks, and the frequency with which they occur, remain unknown. The large 
disk mass suggests the possibility of very large dust fluxes and the commonality of binaries 
suggests that the frequency could be high. Other examples of stars with circumbinary disks 
are HD 44179 (the central star of the Red Rectangle nebula; Jura, Turner, & Balm 1997), 
IRAS 09425-6040 (Molster et al. 2001), and SS Lep, 3 Pup, and BM Gem (Jura, Webb, & 
Kahane 2001). 

If the AMOR point source is (or was) an AGB star, the star would currently be a young, 
hot, high proper motion white dwarf. The proper motion would point back to the vicinity 
of the AMOR source, after allowing for the uncertain aberration of the dust particles. 

6. DUST FROM YOUNG STELLAR OBJECTS 

The accretion of material onto a forming star is accompanied by a bipolar outflow of 
material. These outflows appear as both highly coUimated, supersonic {v ~ 100-200 km s-^) 
jets and as large-angle molecular outflows (with v < 25kms~^). See Reipurth & Bally 
(2001), Konigl & Pudritz (2000), Shu et al. (2000), Eisloffel et al. (2000), and Richer et 
al. (2000) for reviews. 

The mechanism for launching the outflow is not yet well understood. If the material 
originates in the disk and if the density is high enough, then large grains can be entrained in 
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the outflow, assuming there has been sufficient time for growth to large sizes. As the density 
decreases, the grains decouple from the gas and retain their launching speed, which could 
be ~ 10-500 km s^^. Chugai (2001) has appealed to the ejection of large grains from young 
stellar objects (YSOs) to explain the early AMOR detections of interstellar meteors (Taylor, 
Baggaley, & Steel 1996). He finds that the flux supplied from YSOs is a factor of 30 below 
that claimed by the AMOR group. 

Shocks form when the jets run into the ambient ISM and grains could be destroyed in the 
shocked regions (known as Herbig-Haro objects). Observations have yielded contradictory 
results on this point. Beck-Winchatz, Bohm, & Noriega-Crespo (1996) and Bohm & Matt 
(2001) find that the gas-phase Fe abundance in Herbig-Haro objects is generally close to 
solar, indicating very efficient grain destruction. Mouri & Taniguchi (2000), on the other 
hand, find the gas-phase Fe abundance to be fa 20% solar, indicating only a modest amount 
of grain destruction. Grain destruction is not expected for the slower, wide-angle outflows. 
However, these outflows may consist of ambient material swept up by the jets (Konigl & 
Pudritz 2000), in which coagulated grains might not be present. Thus, it is not yet clear 
whether or not YSO outflows contain large grains. 

We can make an optimistic estimate of the dust luminosity by assuming that ~ 1% of 
the outflowing mass is in large grains, as we did for AGB stars. The mass loss rate can be 
as high as ~ IQ-^M^yr-^ (Eisloffel et al. 2000), yielding a flux 



where we have assumed a = 25;um. For ejection in a jet, this is ~ 100 times smaller than 
the optimistic estimate for an AGB star at the same distance. However, the space density of 
YSO's is larger than that of AGB stars; TW Hydrae, a T Tauri star, is only 59 pc from the 
sun, and is associated with dozens of other young stars at similar distances. There is good 
evidence that TW Hydrae has a ~ lOOkms"^ wind Herczeg et al. (4); Wilner, Ho, Kastner, 
& Rodriguez (19). TW Hydrae may be associated with Gould's Belt Makarov &: Fabricius 
(11). While the (optimistic) flux estimate given here is lower by a factor of order 100 than 
is realistically detectable by AMOR. TW Hydrae and stars associated with it might be seen 
by more sensitive future radar systems. 

The outflow may last ~ 10^ yr, 10 times longer than the AGB star lifetime. In the 
following stage of evolution, the YSO disk is dispersed, perhaps by photoevaporation. In 
this phase, the mass loss rate could be as high as ~ 10"'^ Mq yr~^ (Hollenbach, Yorke, & 
Johnstone 2000), yielding a maximum dust luminosity 10 times lower than for the bipolar 
outflows. 




(100) 
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Given the potentially large dust luminosities, YSOs could be likely candidate sources 
for extrasolar meteoroids. However, the young stars listed in Table 4 are probably too old, 
given their proximity. For example, if disk dispersal is complete at an age of 3 Myr, then dust 
from the YSO phase of /?Pic has been traveling for ^ 9 Myr. At a speed of lOkms^^, the 
dust would now be ^ 90 pc from f3 Pic, but the Sun is only ^ 20 pc from j3 Pic. If, however, 
disk dispersal continues until an age of 10 Myr, then (3 Pic may appear as a YSO dust source 
today. 

Some nearby, young clusters could produce large total fluxes distributed over a wide 
area on the sky. Such clusters include TW Hydrac {d ^ 55 pc, t ~ 10 Myr; Zuckcrman & 
Webb 2000), Upper Scorpius {d ~ 145 pc, t ~ 1-10 Myr; Spanglcr ct al. 2001), Chamaclcon 
la and lb {d ~ 140-150 pc, t ~ 1-40 Myr; Spangler et al. 2001), and Taurus {d ~ 140 pc, 
t ~ 10-40 Myr; Spangler et al. 2001). 

7. DISCUSSION 

Our calculations suggest that there arc a number of nearby sources of lOfim or larger 
particles that yield sufficient fluxes to be detectable by ground based radar systems. Here 
we suggest ways to optimize radar systems for meteor detection. 

Meteoroids with a < lO/^m cannot travel through the ISM for appreciable distances. 
Furthermore, it appears that the flux of meteoroids decreases with increasing meteoroid size. 
It follows that if one is interested in detecting particles from an identifiable source, large 
radar collecting area is more crucial than large radar power, once the power is sufficient to 
detect a ^ 10//m particles. 

The returned radar power falls off as for radars that detect coherent emission 

from a substantial fraction of the meteor trail (such as AMOR), or as 1/i?^ for meteor 
head detectors such as Arecibo. The minimum range to the meteor trail is given by the 
height of the trails (of order 100 km), while the maximum range is of order 1000 km due 
to the curvature of the Earth. Increasing the mean range by a factor of 10 will increase 
the collecting area, and hence the flux, by a factor of 100, but requires an increase in radar 
power by a factor of 10^ to lO"', depending on the type of radar employed. While the cost of 
a high power radar transmitter increases rapidly with increasing power, it is clearly helpful 
to maximize the transmitted power. 

Another way to increase the radar range, and hence the collecting area, is to increase 
the gain of the antenna; the received power scales as the product GtGr. However, there are 
limits to the extent to which one can increase the antenna gain. If the radar is like AMOR, 
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the width of the beam must exceed the length of a typical meteor trail, else the number of 
coherently emitting electrons will drop. The angular width of the trail is of order Hp/R, 
which ranges from 0.06 to 0.006 radians, or 3° (in RA, if the beam looks due south or north) 
at the zenith to 0.3° at the horizon. Optimally, the beam will cover the sky from the zenith 
to the horizon, so the antenna gain should be no larger than ~ 1000. 

Selecting a gain of this magnitude matches the beam size to the length of the meteor 
trail, but limits the geometric area Aq of the antenna, partially defeating the purpose, 
which is to maximize the collecting area. This is the case with AMOR, where the choice 
of a narrow beam was dictated by the need for high precision measurement of the meteor 
position, but at the cost of collecting area. To get around this problem, we propose building 
an array of antennas, all powered by the same radio generator. The beams of the array 
would all have widths of order 3°, but the beams would be directed around the points of the 
compass. In principle one could have ~ 100 such beams emanating from a single facility. 
The duration of a radar pulse is typically a microsecond, while the time between pulses is 
of order a millisecond. (The time to travel out and back 1000 km is about 6 milliseconds; in 
that time the meteoroid will travel a distance ~ 240 m.) One could send out pulses, directed 
at different points around the compass, separated by 10 milliseconds. The receiver would be 
turned off when the transmitter is on, about 10~^ of the time. Return pulses from trails at 
different compass points would overlap in time, but given sufficient receiver gains, cross talk 
between the receivers should be small enough to be acceptable. 

If one employs a sufficiently powerful radar, one can live with a low antenna gain and 
still achieve large collecting area. However, the need for high precision determinations of the 
meteor positions remains. This problem can be addressed by using separate transmitting and 
receiving antennae. The transmitter can have a broad beam, with a low gain but covering a 
large fraction of the sky. An antenna array can then be used as an interferometer to locate 
the meteor on the sky. 

Either a multi-beam radar or single broad beam radar with an interferometer could have 
a collecting area Acoi of order 800 km^, if the typical range were i? ~ 150 km. Current distant 
early warning (DEW) radar may be interesting in this context. The existence of tens of such 
systems, with an aggregate collecting area of ~ 10, 000 km^ only adds to their attractiveness 
as extrasolar meteor detectors. 

The collecting area grows dramatically, up to 35, 000 km^, if the radar can detect 10//m 
particles out to the horizon. This is 4000 times the collecting area of the AMOR radar. 
AMOR sees about 10 extrasolar meteors every day; the proposed detector would see up to 
40, 000 extrasolar meteoroids of size 35;um a day, and about 10^ lOyum extrasolar meteoroids 
per day. This suggests that one of the major constraints is handling the data, since there 
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would be of order 10^ interplanetary dust particles detected per day. 

Such a radar would be capable of detecting meteoroids from a handful of nearby debris 
disks; table 4 lists four nearby young Gliese stars that might have fluxes exceeding 6 x 10~^ 
particles per year per square kilometer. This large area radar system would also detect 
meteoroids from (former) AGB stars with fluxes as low as 2 x 10"^^ cm~^ s~^; in the optimistic 
case that a substantial fraction of emitted particles have a > 25/im, sources as distant as a 
kiloparsec could be seen. From table 5 this could be several hundred stars. The number of 
possible YSO sources could be somewhere between these two estimates. 

How feasible is such a radar? We have already mentioned the current DEW systems. 
However, we note that moderate enhancements to an AMOR type radar are all that are 
needed. The product of the transmitter and receiver gain could be enhanced by a factor 
of 5 — 10, while the power could be enhanced by a factor of 20, to 2 megawatts, for a 
total increase in received power of a factor of ~ 200. This would allow for the detection of 
a — 10//m size particles at ranges of ~ 500 km. Further progress could be made by increasing 
the operating wavelength; a factor of two increase, to A = 2000 cm would be sufficient to 
detect 10/im size particles out to the horizon. 

The flux of interstellar meteoroids of different mass appears to follow a power law 
mfm oc m~^'^, over five decades in mass. This may be compared to the Dohnanyi law, 
fm oc m~", with a — 11/6 1.83; in the case of interstellar meteoroids, a ^ 2.1. We 
note that the observed size distribution of asteroids has a 2 for objects larger than 
a — 2.5km, and a ^ 1.4 for a between 2.5 and 0.1km (Ivezic et al. 2001). However, the 
binding energy per gram can be expected to vary with size, a fact wc appealed to in explaining 
the observation that meteoroids of very difference sizes ablate at the same height in the 
Earth's atmosphere. In our view, the rough agreement between the values of a might reflect 
a common origin in a collisional cascade, but is more likely a coincidence. 

Krivova & Solanki (9) recently suggested that interactions with a Jupiter mass planet 
could eject dust grains with velocities of order 10 — 70kms~^. They cited this result as 
support for the identification of the point source found by Baggaley (2000) with j3 Pic. Our 
analytic calculations show that only a tiny fraction of ejected particles will have such high 
velocities, a result supported by direct numerical integrations. 

8. CONCLUSIONS 

We have provided simple analytic estimates for the minimum size of radar detected 
meteoroids as a function of radar power and antenna gain. We give very rough estimates of 
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the collecting area of radar detectors, and hence of fluxes of interstellar meteoroids, assuming 
that the claimed detection rates are correct. The fluxes of satellite and radar detected 
interstellar meteoroids appear to lie along a power law, mfm oc m~^'^, corresponding to a 
a — 2.1, or a size scaling 7 = 4.3, in contrast to the Dohnanyi value 7 = 3.5. 

We examine three possible sources of large, a > 10//m size extrasolar meteoroids. In 
descending order of dust luminosity, they are AGB stars, young stellar objects (with winds 
and/or jets), and debris disks. (Supernovae can supply similar mass fluxes of particles, but 
we expect the particles to be smaller.) We show that such large particles can travel in 
straight lines for tens of parsecs through the interstellar medium, which in principle will 
allow their sources to be traced. 

Since there are ~ 2 AGB stars within lOOpc of the sun at any time, we expect that there 
will be ~ 2 sources of extrasolar meteoroids, with fluxes at Earth approaching ~ 8 yr~^ km~^, 
on the sky at all times. Such a star might be responsible for the "point source" reported 
in Baggaley (2000). The AMOR radar, according to our estimates, is marginally capable of 
detecting such a small flux. The star itself would currently be a hot white dwarf with a high 
proper motion, within ~ 100 pc of the sun. 

There is one known T Tauri star within 100 pc of the sun, TW Hydrae. This star has 
a high velocity outflow (~ lOOkms"^) from a dusty disk. It is very likely expelling 10 to 
100/xm size particles; unfortunately, the flux is clearly too small, by a factor of ~ 100, to be 
detected by current radar systems. 

Similarly, currently known debris disk systems provide fluxes that arc too small to be 
detected at present. We have shown that the system imphcated by Baggaley (2000), /3 Pic, 
is unlikely to provide a particle flux as large as that associated with his point source. We 
also show that the location of the apparent source on the sky is inconsistent with particles 
ejected from /3 Pic by gravitational interaction with a Jupiter mass planet in the system; we 
showed that typical ejection velocities in that case are of order 1 kms""*^. Even if the ejection 
velocity is assumed to be 30kms~^, the most favorable value, the apparent position of the 
source does not match the observed position (see Fig. 8). Since the particle flux, ejection 
velocity, and position on the sky do not match the observations, we conclude that the point 
source is not associated with gravitational ejection from the debris disk of P Pic. 

However, we would like to stress that all three types of sources, AGB stars, YSOs, and 
debris disks, should be detectable by future ground based radar systems. The pioneering 
AMOR system may have already detected meteoroids from an AGB star. Modest improve- 
ments in this type of ground based radar system, in particular better sky coverage, should 
allow for the detection of multiple examples of each type of source. An alternate possibility 
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is to piggyback on already existing radar employed in distant early warning systems. 
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bourg, France, and of NASA's Astrophysics Data System. 



A. APPENDIX 

In this appendix we rederive some of the relations used in section 2, employing the 
classical theory of meteors. In the meteoroid frame, the air flow has a kinetic energy flux 

F = Ipav'. (Al) 

Following the meteoritics hterature, we define a shape factor A to be the ratio of the area 
Am of the meteoroid divided by the two-thirds power of the volume (m/p), or 

^= (m/p")V3- ^^2) 

For a sphere, A — (Qtt/IG)^/^ 1.2. The kinetic luminosity seen by the meteoroid is then 
AjnF — ^{m/py^^paV^; assume that a fraction A of this kinetic luminosity goes toward 
ablating the meteor. Then the time rate of change of the binding energy is 

dE KAfmV'^ 3 



Typical estimates for A are around 0.5. It follows that the rate of ablation is 

dm KA /m^^^^ 



. . , PaV, (A4) 
dt 2C \p J ^ ' 

where, the reader will recall, the binding energy E is related to the meteoritic mass m by the 
heat of ablation (, E = C,m. From equation (8), the number of ions produced per centimeter 
along the path is 

9 = (A5) 

vpL dt 

where is the mean molecular weight of the meteoroid (recall that P is the number of ions 
produced by each meteoroid atom). Combining this with equation (A4), we find 

^KA (mV'^ 2 
? = - Pav'. (A6) 
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This is known as the ionization equation. 

The maximum hne density along the track can be found by differentiating the ionization 
equation with respect to time: 

Idq 1 dp{v) 2 dm 1 dpa 2 dv , , 

q dt P dt 3m dt pa dt v dt 

Note that f3{v) ~ f", so that both the first and last terms on the right hand side are 
proportional to 

(A8) 

vdt ^ ^ 

We assume for the moment that this is smaller than either the mass or density variations. 
Then the maximum value of q occurs when 

2 dm 1 dpa 
3m dt Pa dt 

or 



(A9) 



where Hp is the density scale height. Thus 

W = (All) 

In words, the maximum line density is given by spreading the meteoroid over approximately 
a scale height and accounting for the number (3 of atoms that are ionized for each atom in 
the meteor. 

Next we justify the neglect of terms proportional to the derivative of the velocity. The 
momentum equation for the meteoroid is 

dv ( mi \ 

m— ^mg- CdA \ — \ pav'', (A12) 

where Cd is the drag coefficient. The mean free path at meteor heights is of order 10 cm, 
much larger than the typical radius of the meteoroids we are interested in. In that case 
Cd = 1. The velocity derivative becomes 

Idv g fmV^^ 

-— = - -AI — ] pav m A13 
V dt V \ P J 

We note that Hp c^/g where c is the sound speed; then 

— fti - A14 
H„ \c/ V V 
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so that the first term on the right hand side of (A13) is neghgible compared to v/Hp\ recall 
that the meteoroids are highly supersonic. Next we compare the velocity gradient in (A7) 
(including only the second term on the rhs of (A13)), 

—^A[-] t^r^a^p^v pa? A15 

V at \ p J m p a Hp 

2 dm AA (mV'^ pav^ 

= 1-1 — (A16) 



to 



The ratio is 



3m dt 2( \ p J 

12 c 3 / c \ r V \2 /lo 



A 4 V 10" cm^s-^ y V40kms-V V A 

where we took n = 2. 



(A17) 



McKinley gives estimates for A that arc of order unity. For these values, the velocity 
derivative is much smaller than the mass derivative. We conclude that the logarithmic mass 
derivative is comparable to or larger than the logarithmic velocity derivative. 
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Fig. 1. — Measured cumulative fluxes of extrasolar meteoroids. The data are from the Ulysses 
and Galileo satellites (crosses joined by a solid line), the Arecibo radar (solid squares) and 
the AMOR radar (solid triangle). The open triangle represent the "point source" seen by the 
AMOR radar. The dashed line represents the mfm oc m~^-^ scahng first noted by Landgraf 
et al. (2000) and given in equation (5). The upper axis gives the particle radius a in microns, 
assuming that the meteoroids have a density of 3gcm~^. 
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Fig. 2. — Measured differential mass fluxes, m ■ mfm, with units of gcm~^s~^, of extrasolar 
micrometeoroids. The long-dashed line is the Landgraf et al. (2000) scaling, while the 
dotted line is the limit found by assuming that half of all the metals in the local ISM (with 
n = 0.1 cm~^) are in the form of micrometeoroids. The solid squares are the fluxes reported 
by the Arecibo radar. The filled triangle is the fiux estimated from the spatially distributed 
fiux in Fig. 2 of Baggaley (2000). The open triangle is the fiux in Baggaley's "point source" . 
The horizontal dotted line is the upper limit on the fiux given by assuming that half of 
all the heavy elements are in dust grains. The dashed line joining open circles is the fiux 
corresponding to the ISM dust size distribution found by Kim and Martin (1995), normalized 
to the dotted line. This figure suggests that the bulk of the mass in dust consists of particles 
with masses of order 10~^^ g, or sizes of order two tenths of a micron. 



-53- 




Fig. 3. — The drag distance i^drag (short dash), gyroradius (sohd), and grain destruction 
distance D^est (long dash) versus grain speed v for four grain sizes, as indicated, and assuming 
supersonic grain speeds and p = 3.5gcm~^, riu = lcm~^, U — 0.5 V, and B = b/iG. 
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50 100 150 200 

d (pc) 

Fig. 4. — The upper and lower current ages i of a Vega-like star for which we would observe 
grains emitted when the star's age was between tcv and 2tcr, versus the distance d to the 
star. 
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Fig. 5. — The geometry of hyperbolic motion in the two body problem. The scattering 
center (the planet) is at the origin; the test particle travels along the curved line, with an 
initial impact parameter given by s. Its closest approach to the planet occurs as it crosses 
the apsidal line (choosen here to be the x-axis), when it is a distance q from the planet. 
The ingoing and outgoing asymptotes are denoted by dotted hnes. The angle between the 
ingoing and outgoing asymptotes is 7. 
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Fig. 6. — Ejection velocities of small particles from a system with a solar mass star and a 
Jupiter mass planet at 5.2AU. 



-57- 



90 



60 - 



30 - 



-30 - 



-60 - 



-90 



I I I — TT 



I I I — I I I I I 



I I I I I 



I I I I I — I I I 



□ 



J I I I I I I I 1 I I I I I I I I I I I I I I I I L 







360 



300 



240 



180 
A (deg) 



120 



60 







Fig. 7. — Triangles: apparent locations of the dust flux from the Gliese stars in Table 4. 
Open box: Direction of the solar apex. Coordinates are ecliptic longitude (A) and latitude 
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Fig. 8. — The box marks the central location of the discrete source observed by Baggaley 
(2000) and the dashed line indicates its extent. Triangles mark the apparent location of a 
dust stream from /5Pic for values of the ejection speed v^] ranging from 2 to 40kms~^ in 
steps of 2kms~^. For a few cases, Vej and the speed of the dust with respect to the Sun, v^^q 
are indicated. The ecliptic coordinates of /3Pic itself are (A, /3) = (82.5°, —74.4°). 
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Table 1. Radar Parameters 





Arecibo 


AMOR 


Pt 


2MW 


lOOkW 


Pn 


10"^ erg/s 


1.6 X 10"'^ erg/s 


Gt 


10^ 


430 


Gr 


10^ 


130 


Ro 


100 km 


230 km 


A 


70 cm 


1145 cm 


Beam dimensions 


1 ° 1 ° 

6^6 


3° X 18° 


Ag 


0.9 W 


50 km^ 


A-col 


10-2 W 


Skm^ 


meteor properties 


Qmin 


lO^cm-^ 


7 X 10^ cm-^ 




3 X lO-'^ gm 


2 X 10-^ gm 




6 lira 


25 /im 
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Table 2. Idealized Interstellar Environments 



Phase 


T 

-'gas 








K 






CNM 


100 


30 


1.5 X 10-3 


WNM 


6000 


0.3 


0.1 


LB 


10^ 


5 X 10-3 


1 



Table 3. Threshold Distances for Observing Vega-Like Stars 



A,:ol 






Oct 


km^ 


Myr 


pc 


//m 


10^ 


3 


26 


17 


10^ 


10 


8.6 


10 


10^ 


3 


108 


35 


10^ 


10 


36 


20 
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Table 4. Gliese Star Dust Fluxes 



Gl# 


Name 


d 










^^,0 








(3' 






pc 


Myr 


kms~^ 


kms~^ 


kms~^ 


kms~^ 




yr~^ km~^ 


deg 


deg 


68.0 


DM+19 279 


7.5 


69.8 


34.5 


Ct A IT 

-24.7 


r\ A 

2.4 


Ad r" 

42.5 


A A 


3.3E-5 


62 


33 


71.0 


rCet 


3.6 


170 


18.6 


29.4 


12.7 


37.0 


3.3 


6.3E-6 


94 


-54 


111.0 


t'- En 


14.0 


39.6 


-26.2 


-16.9 


-13.0 


33.8 


6.8 


1.5E-4 


262 


33 


121.0 


En 


26.4 


215 


19.4 


9.96 


-1.09 


21.9 


11.6 


3.4E-8 


101 


-18 


144.0 


eEri 


3.22 


59.7 


-3.01 


7.22 


-20.0 


21.5 


4.1 


l.lE-4 


184 


8 


167.1 


7 Dor 


20.3 


107 


-22.9 


-18.1 


-13.8 


32.3 


8.4 


1.2E-6 


261 


38 


217.1 


CLep 


21.5 


69.8 


1/1/1 

-14.4 


-11.2 


-8.45 


20.1 


10.9 


3.7E-6 


260 


38 


219.0 


p PlC 


19.3 


12.0 


-10.8 


-16.0 


-9.11 


21.4 


10.0 


1.5E-3 


268 


51 


248.0 


a ric 


30.3 


194 


-35.4 


-19.3 


-10.1 


41.6 


9.0 


l.OE-7 


268 


25 


297.1 


B Car 


21.4 


O 1 A 

31.4 


19.8 


-15.1 


-34.0 


Ana 

42.2 


7.5 


l.lE-3 


140 


A C\ 

48 


321.3 


5Vel 


24.4 


351 


11.6 


-1.15 


-5.00 


12.7 


14.6 


2.3E-9 


107 


17 


364.0 


DM-23 8646 


14.9 


69.3 


-39.9 


-25.4 


6.27 


47.7 


5.9 


1.4E-5 


290 


17 


448.0 


;5Leo 


11.1 


135 


-20.1 


-16.1 


-7.80 


26.9 


6.8 


1.2E-6 


268 


36 


557.0 


<T Boo 


15.5 


123 


2.06 


15.9 


-5.17 


16.9 


10.1 


5.7E-7 


165 


-38 


673.1 


DM-24 13337 


25.7 


34.5 


-35.8 


-13.1 


-11.6 


39.9 


8.5 


6.4E-4 


261 


20 


691.0 


/i Ara 


15.3 


58.8 


-13.7 


-8.39 


-3.97 


16.5 


10.1 


9.7E-6 


270 


26 


721.0 


a Lyr 


7.76 


152 


-16.1 


-6.33 


-7.76 


18.9 


6.7 


l.lE-6 


256 


26 


820.0 


61 Cyg 


3.48 


48.0 


-93.8 


-53.4 


-8.33 


108 


1.9 


l.OE-3 


278 


23 


822.0 


5Eql 


18.5 


80.0 


5.72 


-28.6 


-10.3 


30.9 


8.2 


4.3E-6 


18 


78 


881.0 


aPsA 


7.69 


83.1 


-5.70 


-8.22 


-11.0 


14.9 


7.6 


5.7E-6 


235 


54 



*Age of star from eq. (96) 

'^Space velocity of the star relative to the Sun; U is positive towards the Galactic center, V is positive 
along the direction of the Galactic rotation, and W is positive towards the North Galactic Pole. 

■^Dust flux calculated using equation (92) with a limiting grain size a = 10/xm, a flat distribution of 
ejection velocities between 0.5kms~^ and 3kms~^, Tcr = 3 Myr, and /survive = /beam = /ej = 1- 

'^Apparent location of dust stream on the sky; A is ecliptic longitude and /? is ecliptic latitude. In most 
= 1 kms ^ is adopted; for Gl 219.0 (/SPic) Vej — 3 kms ^ is adopted, since the dust will not yet 
have reached us if Vei = 1 kms~^. 
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Table 5. AGB Star Dust Fluxes 











9' 


pc 










100 


7.99 


1.84 


24.3 


2.50E-4 


200 


0.706 


14.7 


34.3 


2.83E-3 


500 


2.86E-2 


230 


54.2 


6.99E-2 



"^distance to AGB star 
''flux from eq. (99) 

'^number of AGB stars within distance 
d from Sun 

'^threshold grain size for traceability to 

source 

^fraction of emitted grains that must 
have a > at in order for flux at Earth 
to equal 2 x 10~^yr~^km~^ 



